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1. INTRODUCTION 

The work supported under NSG-6022 entitled "Rocket Study of Auroral 
Processes" commenced in 1978 and terminated June 30, 1981. Under NSG-6022 
the University of New Hampshire built particle experiments for three rocket 
flights. These were a joint University of New Hampshire - University of 
Minnesota Terrier Malemute Flight, NASA 29:014, launched January 27, 1980; 
University of Minnesota, Echo V Flight, launched November 13, 1979; and the 
Norwegian Corobier Ferdinand 50 Flight, launched January 24, 1980. This 
grant also supported the analysis of data obtained from earlier rocket 
flights, Echo III, Flight 18:165, Flight 18:1004 and Flight 18:1005, and 
theoretical studies. 

The form of this report will discuss each of the items in the work 
statement in detail starting first with the data analysis of earlie ' flights 
and theoretical studies. 

2. ANALYSIS OF DATA - EARLIER FLIGHTS 


2,1 Auroral Electrodynmics 

2.1.0, Theoretical paper - "Modulation of Terrestrial Ion Escape 
Flux Composition (by Low-Altitude Acceleration and Charge 
Exchange Chemistry)", T.E. Moore, J. Geophys. Res., 8_5 
2011, 1980. 

Abstract 

Motivated by recrnt observations of highly variable hot 
plasma composition in the magnetosphere, control of the 
ionospheric escape flux composition by low-altitude 
particle dynamics and ion chemistry has been investiga- 
ted for an e~, H"**, 0 + , ionosphere. Theoretical polar 
wind results of Banks and Holzer and observations of 
highly nonthermal ion escape have been used as guides. 

It is found that the fraction of the steady state es- 
cape flux which is 0"*" can be controlled very sensitively 
By the occurrence of parallel or transverse ion accel- 
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eratlon at altitudes bt'low the altitude where the 
neutral oxygen density falls rapidly below the neutral 
hydrogen density and the ionospheric source of 0 + tends 
to be rapidly converted by charge exchange to 11+ . The 
acceleration is required both to overcome the gravita- 
tional confinement of 0+ and to violate charge exchange 
equilibrium so that the neutral hydrogen atmosphere 
appears 'optically' thin to escaping 0+. Constraints 
are placed on the acceleration processes, and it is 
shown that 0+ escape is facilitated by observed iono- 
spheric responses to magnetic activity, 

2.1,1 Paper published using Plight 18:165 data - "Rocket Observa- 
tions at the Northern Edge of the Eastward Electrojet" , 

L.J. Cahill, Jr., R.L. Arnoldy, W.W.L. Taylor, J. Geophvs. 
Res., 85, 3407, 1980. 

Abstract 


A Nike-Tomahawk rocket was launched north over two quiet 
late evening arcs in March 1975. A northward magnetic 
disturbance x«is observed on the ground under the rocket 
trajectory. South of the arcs the northward electric 
field was 60 mV/m, indicating strong westward plasma 
flow. An eastward electrojet current layer was pene- 
trated in the upward flight. Precipitating electrons 
were observed over each arc. The electric field de- 
creased to below 20 mV/m over each arc and recovered to 
50 mV/m between the arcs. Using the observed electron 
flux and a model of the ionosphere, the Hall and Peder- 
sen conductivities were calculated. These conductivi- 
ties were used, with the observed electric field, to 
calculate the eastward and northward components of the 
horizontal ionospheric currents. The eastward current 
calculated south of the first arc agreed well with the 
observed eastward electronjet current. The power dlsi- 
pated by the Pedersen current S p | E | 2 was also calculated 
and compared with the power carried by the precipitating 
electrons. The Joule power decrease abruptly in the 
auroral arcs, as the precipitating electron power ln~ 
creased. The total dissipated power was approximately 
th'e same inside the arcs, between them, and south of 
the luminosity. North of the aurora the electric field 
and dissipated power remained low. Field-aligned cur- 
rents carried by the observed electrons were about a 
factor of 3 lower than those inferred from the magnetic 
field measurements. Likewise, current continuity argu- 
ments to keep the auroral current system divergence 
free required a larger field-aligned current than that 
obtained from the particle detectors. 

2.1.2 Paper published using Flight 18:1005 data - "Ionospheric 
Electrical Currents in the Late Evening Plasma Research", 
D.A. Behm, F. Primdahl, L.J. Zanetti, Jr., R.L. Arnoldy, 
and L.J. Cahill, Jr., J. Geophys. Res., 84, 5339, 1979. 

j 
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than that from measured precipitating electrons. There were 
also discrepancies between the Ion drift measurement and the 
other measurements. We conclude that the various techniques 
are subject to large error o and that measurements of this im- 
portant quantity be used with caution. 

2.2 Electron Gun Dynamics 

2.2.1 Paper published using Echo ill data - "The Hot Plasma Environ- 
ment and Floating Potentials of an Electron-Beam-Emitting 
Rocket in the Ionosphere", R.L. Arnoldy and J.R. Winckler, J. 
Geophys. Res., 86 , 575, 1981. 

Abstract 


With an extensive array of particle sensors the plasma environ- 
ment surrounding the Echo III accelerator payload is studied. 

From measurements of the thermal ion spectrum, negative payload 
potentials referenced to the unpertured ionospheric plasma are 
obtained. Multiple detectors determined the electron popula- 
tion from a fraction of an eV up to AO keV energy. An intense 
electron population extending up to the energy of injected elec- 
trons by the accelerator is produced when the acceleration is 
turned on. The energetic tail of this population is returned to 
the payload principally from directions in which the beam was 
fired, with the most intense fluxes coming from the atmosphere 
for downward injections of the beam. The atmospheric-scattered 
beam and secondary electrons are called 'quick echoes'. Elec- 
trons of energy less than several lceV down to the detector thres- 
hold (0.1 eV) are called the suprathermal component. These elec- 
trons are produced isotropically around the payload during gun 
firings and decay away in approximately 32 ms. The largest di- 
rectional intensities of this component are observed at the 
higher altitudes. Quick echo electrons are also observed to pro- 
duce suprathermal electrons when they encounter the payload. 

The mechanism by which the suprathermal electrons are produced 
is discussed but remains unknown at the present time. The hot 
electrons surrounding the accelerator payload during gun injec- 
tions bring sufficient charge to the payload to neutralize it 
provided the loss of charge by secondary production on the pay- 
load skin is small, presumably owing to a positive payload float- 
ing potential during injection. Since the hot population exists 
for tens of milliseconds after the gun turn off, it results in 
driving the payload up to 4 volts negative during this time. 
Secondary production on the payload skin apparently prevents 
much larger negative potentials after gun turn off when the 
payload is immersed in the hot population. Quick echo electrons 
creating suprathermal electrons around the payload also drive 
the paylaod to a few volts negative. This is an important con- 
sideration when discussing the origin of the suprathermal elec- 
trons, since electrical discharge and beam plasma discharge 
mechanisms do not apply to the quick echo beams. 


3. TERRIER-MALEMUTE FLIGHT- NAJA JL9:014 


3.1 Description of Payload and Flight 


>pr~: 
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This rocket payload was launched up to an altitude of ^ 400 km over a quiet, 
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pre-midnight, auroral arc. It was Cully instrumented to measure AC/DC 
fields, plasma, and energetic particles to Investigate the electrodynamics 
of the aurora and its coupling to the magnetosphere. In addition to these 
passive measurements the flight contained an active experiment to study 
the plasma dynamics of firing 80 milliamperes of 60 e.V argon ions into 
the ionosphere under auroral conditions. 

AC/DC electric and magnetic fields were measured by the University of 
Minnesota. In actuition the. University of Minnesota built the argon 
plasma generator. Bernard llaQsler of the Max Planck Institute, Catching, 
provided a search coll magnetometer. Dr. Edward Szuszczwicz of the Naval 
Research laboratory diagnosed the electron plasma with a pulsed plasma 
probe. The University of New Hampshire made the ion plasma measurements, 
topside photometer observations, and the energetic ion and electron mea- 
surements. 

Multiple electrostatic analyzers measured the electron and ion distribu- 
tion function from a few eV to 30 keV every 200 ms in 64 energy steps. 

In addition, with a peak tracking analyzer, 0.8 ms time resolution was 
obtained at the energy of the peak in the electron spectrum, 

Finally, the University of New Hampshire ion plasma experiment measured 
the thermal ion population with four detectors positioned on the spinning 
payload to obtain a distribution function approximately once per second. 

The detectors were small geometric factor, electrostatic analyzers, sweep- 
ing the ion spectrum from 0 5 eV in 32 steps once every 25.6 milleseconds. 

3.2 Flight Summary 

The flight was highly successful. All experiments functioned properly 
with the exception of a great .deal of interference in the search coil 
data. Vehicle performance was nominal, however, the electronic’ payload 
timer was late by about 34 seconds. This delay in function timing (see 
Wallops Report "Post Flight Investigation of the Terrier-Malemute 29:014 
UE Delayed Mission Events, May 1980) resulted in 1) loss of 100 km up- 
leg electrojet data, and 2) failure to recover payload. 

The flight was in a diffuse aurora on the upleg and, according to ground 
data, over the eastward electrojet. Just beyond apogee the flight en- 
tered discrete aurora and crossed the Harang discontinuity into the west- 
ward electrojet. On the downleg of the trajectory and during the period 
of ion plasma release the flight was in transient discrete aurora. 

3.3 Results 

3.3.1 Paper presented at Fall A.G.U. 1980 - "Preliminary Results of 
Plasma Gun Operation in the Auroral Ionosphere", L.J. Cahill, 

Jr., R.L. Arnoldy , and D. A. Behm, EOS, 61, 1082, 1980. 

Abstract 

A sounding rocket instrumented with a full set of particles and 
fields detectors was launched from Poker Flats into the evening 
auroral ionosphere in January 1980. A 60 volt Argon ion gun on 
the rocket was operated intermittently during the last third of 
the flight. Immediately when the ion gun was turned on a number 
of effects were noted in the detectors. The signal from the 
electric field booms was saturated during part of each spin 
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period while the gun was on but displayed the normal sine 
wave modulation during the remainder of the period. A strong 
2 to 4 kHz signal was received In an AGO receiver connected 
to the set of electric field booms. This signal was also 
spin modulated with large amplitude during only part of each 
spin period and amplitude lower by at least one order of mag- 
nitude during the remainder of the' period. There was also a 
large Increase In the waves observed with a 1 to 10 mHz swept 
frequency receiver. 

The gun operation also produced dramatic changes in the charged 
particle fluxes. Two new populations of electrons appeared. 

The first was an enhancement betwee 100 and 1000 eV with the 
peak, flux at 300 eV. The flux increased from 10 8 to 10 9 (cm 2 
sec sr keV)" 1 at 300 eV. This population was Isotropic with 
as many electrons coming up the field line as down. The second 
population was below 150 eV, increasing by a factor of five for 
electrons streaming down the field line and decreasing for parti- 
cles coming back up. 

3.3.2 Invited paper presented at "Chapman Conferences on the Formation 
of Auroral Ares", 11,1. Amoldy, College, Alaska, July, 1981. 

Abstract 


Auroral electron and ion precipitation measured by sounding 
rockets and low altitude’ satellites will be reviewed. The two 
observing platforms provide a different perspective of the aurora. 

An attempt will be made to relate the two data sets to establish 
the precipitation signature of an auroral arc. Energy spectra 
and pitch angle distributions of precipitating particles have 
provided the basis for most auroral acceleration mechanisms. 

The altitude and particle mass dependence of these parameters 
as well as their spatial relationship to the visual aurora will 
be reviewed. 

Although optical data have shown the fine structure (temporal 
and spatial) of an arc, improved response of particle detectors 
is just Beginning to reveal this structure in particle precipita- 
tion. For example, intense fluxes of field aligned electrons 
are often observed having a duration on the order of a second as 
seen by rocket detectors. Similar field aligned bursts fluctuate 
at a few Hertz in satellite data. The intent of this paper is to 
summarize the observations that may aid in the formation of an 
auroral arc theory and which ultimately must be explained by such 
a theory. 

This paper is to be published in the "Proceedings of Chapman Con- 
ference on the Formation of Auroral Arcs", A.G.U. Publication, 
Washington, D.C., S.-I. Akasofu, Editor, 1981. (This paper is 
presented in its entirety in Appendix A). 

3.4 Work In Progress 

Much of the data from the flight are still being analyzed, particularly, 
the aspect of combining the AC/DC field, plasma, and energetic particle 
results in a consistent model. Interesting observations we are currently 
pursuing at UNH are the following: 



3' 4.1. Natural Aurora Measurements 

Using large geometric factor detectors capable of obtaining 
good statistics in measuring the weak diffuse aurora parti- 
cle precipitation we find: 

Ions < 1 keV undergo significant heating in the 
direction transverse to B during diffuse aurora. 

Diffuse aurora conics? 

The diffuse auroral electron background shows 
little evidence of acceleration when the dis- 
crete aurora peaked spectum "grows" out of the 
diffuse precipitation. 

Energetic (> 10 keV) ions and electrons having 
a trapped pitch angle distribution show nearly 
an order of magnitude increase in directional 
Intensity when the rocket moves .from the dif- 
fuse to discrete aurora. 

These results suggest AC acceleration mechanisms might be 
important. In addition , the third result above casts some 
doubt on using the trapped pitch angle distribution as evi- 
dence of closed field lines. Why should trapped radiation 
increase on auroral field lines? 

3.4.2 The Argon Plasma Releases: 

Created a down-streaming electron population 
(< 150 eV) much like that often observed at the 
edges of auroral arcs. 

Created a suprathermal population of electrons 
between 100 and 1000 eV energy with a peak flux 
at ^ 300 eV which was isotropic in direction. 

Triggered the release of a short burst of > 15 
keV ions apparently from the trapped radiation. 

Given the fact that the plasma generator vas a very low 
power instrument, it seems imperative that the plasma re- 
leases triggered instabilities in the natural auroral en- 
vironment to produce the effects observed. Finally, the 
results suggest that local ionospheric conditions could 
play an important role in discrete arc formation. 

ECHO V (NASA 12:037) 

4.1 Description of Payload and Flight 

The Echo V payload launched by a Strypi II vehicle out of the Poker Flat, 
Alaska range carried an electron gun capable of injecting 40 keV electron 
bursts into the magnetosphere at a current of about 0.5 ampere. Earlier 
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flights in this program successfully detected electrons that have echoed 
from the conjugate hemisphere. The primary objective of the Echo V flight 
was to employ this technique to probe the magnetospheric configuration 
of the auroral oval. Are the auroral field lines closed, how far into 
the tall do they extend, are they equipotentials? Additional objectives 
were to continue to study the plasma dynamics associated with the gun 
injections. 

The Principal Investigator for this program was Professor John Winckler 
at the University of Minnesota. The University of Minnesota furnished 
the electron gun, neutral gas release experiment and various sensors to 
diagnose the beam and its echo, these included scintillators for x-ray 
and energetic particle detection and photometers. Professor Paul Kellogg 
at the University of Minnesota performed plasma diagnostic measurements 
aboard the spacecraft with a Langmuir probe and wave measurements from 
the ground. Both the University of Minnesota and Professor Neil Davis 
from the University of Alaska conducted optical measurements for beam 
detection from the ground. Finally, the University of New Hampshire 
provided energetic ion and electron detectors and thermal ion sensors. 

Four electrostatic analyzers measuring ions and electrons in the energy 
ranges, few eV to 2 keV and 50 eV to 30 keV, comprised the University 
of New Hampshire's energetic particle package. These detectors swept 
their respective energy range once every 200 milliseconds sampling 64 
energy levels per sweep. The objective of these detectors was to mea- 
sure the auroral precipitation and any perturbation of it due to the 
gun firing. Scattered beam electrons and secondary electrons produced 
by the beam and/or by the beam plasma discharge mechanism have also been 
studied by these sensors. 

The UNH thermal ion sensors were identical to those discussed above for 
the Terriar-Malenmte flight except that they were rotated by a motor 
driven platform to obtain the distribution function since the payload 
did not spin. The motion of gun electrons in traversing the magneto- 
sphere is sensitive to the convection electric fields they experience. 

The prime objective of the thermal ion sensors was to obtain a measure 
of these fields in the ionosphere. 

4.2 Flight Summary 

The flight was successfully launched with its trajectory over a weak 
aurora and just south of a bright auroral band. The UNH detectors 
all collected data, however, just before apogee the energetic parti- 
cle detectors did not sweep the spectrum but were locked at one energy. 
Although this was a "failure" mode, valuable high time resolution data 
were collected as a result of it, No return Echo electrons were mea- 
sured optically by the University of Alaska-Universlty of Minnesota 
experiments. It is believed that the beam electrons were scattered to 
trapped orbits as a result of wave activity on the auroral field lines. 

4.3 Results 

Seme preliminary results on the cross B dimension of the space charge 
region created by the electron gun were presented in an invited review 
paper authored by J.R, Winckler at the"NAT0 Advanced Research Institute 
on Artificial Particle Beams Ultilized in Space Plasma Studies" held in 
Norway, April, 1981. 
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The size of che electrostatically perturbed region transverse to the 
mugnetic field around an electron-beam-emitting payload may exceed a 
kilometer as measured by the reduction in the thermal ion population 
reaching the payload just after beam emission. 

4.4 Work In Progress 

A major area of current study Involves the electron-beam interactions 
with the ionosphere between 100 - 200 km. 

Beams are a strong source of secondary electrons < 1 keV 
in energy (suprathermal population) . The intensity of 
these electrons increas es with altitude. This beam-plasma 
interaction is quenched by the release of neutral gas. 
Fluctuations in the suprathermal electron component can be 
excited at ambient heavy ion gyrofrequencies by the beam. 

The scattering of beam electrons between 100 and 200 km is 
due primarily to beam-neutral collisions. A weak, neutral 
density independent, mechanism is observed above 200 km al- 
titude. 

5 . C0R0B1ER (F-50) 


5.1 Description of Payload and Flight 

Corobier is an acronum for "Coordinated Rocket and jBackscatter Investiga- 
tion of the E-Region^. The program was a Norwegian one with Dr. Jan Holtet 
of the University of Oslo, the Principal Investigator., and Dr. Ray Green- 
wald of the Max-Planck Insitute, Lindau, a Co-Investigator. The general 
objective of the program was to study plasma Instabilities and irregulari- 
ties in the auroral E-region with a combination of ground-based probing 
of the ionosphere with radar and in-situ measurements. Enhanced iono- 
spheric electric fields excite plasma instabilities in the auroral zone 
E-region. Although these waves are thought to be produced by the two 
stream and gradient drift instabilities there are a number of questions 
that remain. 

A host of AC/DC wave experiments and plasma particle sensors were flown. 
Investigators included Holtet, Troym, Soraas and Staden.es from Norway; 
Falthamar and Bostrom from Sweden, Primdahl and Unstrup from Denmark and 
Evans and Arnoldy from the United States. The flight was launched out of 
the Norway Andoya rocket range on January 24, 1980. The UNH experiment 
consisted of thermal ion detectors identical to those that were flown 
aboard Flight 29:014. The objectives were to measure the ion temperatures 
and the E-region convective electric field in aspect of the specific pro- 
gram objectives. 

5.2 Flight Summary 

The UNH detectors worked normally during the flight although they did not 
fully deploy due to the anomalously low spin rate. The low spin rate re- 
sulted in all of the payload doors (except the UNH door) not to deploy, 
hence the flight was a failure. No attempt has heen made to analyze the 
UNH data primarily because no other (except the DC magnetometer) experi- 
ment operated properly and also no support data are available. 
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REVIEW OF AURORAL PARTICLE PRECIPITATION 


R. L. Arnoldy 

Space Science Center* University of New Hampshire, Durham, N. H. 03824 


Abstract . Auroral electron and ion precipita- 
tion measured by sounding rockets and satellites 
will be reviewed. The two observing platforms 
provide a different perspective of the aurora. 

An attempt will be made to relate the two data 
sets to establish the precipitation signature 
of an auroral arc. Energy spectra and pitch an- 
gle distributions of precipitating particles 
have provided the basis for most auroral accel- 
eration mechanisms. The altitude and particle 
moss dependence of these parameters, as well as 
their spatial relationship to the visual aurora 
will be reviewed. Although optical data have 
shown the fine structure (temporal and spatial) 
of an arc, improved response of particle detec- 
tors is just beginning to reveal this structure 
in particle precipitation. For example, intense 
fluxes of field aligned electrons are often ob- 
served having a duration on the order of a sec- 
ond as seen by rocket detectors. Similar field- 
aligned bursts fluctuate at a few Hertz in sat- 
ellite dataT The intent of this paper is to 
summarize the observations that may aid in the 
formulation of an auroral arc theory. 

Introduction 

Particle precipitation associated with dis- 
crete auroral arcs will be reviewed in this pa- . 
per, emphasizing those observations which are 
pertinent to theories of auroral arc formation. 
Because of the large quantity of auroral parti- 
cle precipitation data available, only observa- 
tions by rockets and satellites generally made 
in the evening sector will be considered. Re- 
cent reviews of hot magnetospheric ions and 
electrons associated with optical auroras have 
been made by Hultqvist (1979, 1974). 

Tills review will be organized according to 
observations made in three regions of space: 
Region l accessible to rockets and low altitude 
satellites (few hundred to a few thousand kilo- 
meters) ; Region 2 extending from 4000 *> 8000 Ian 
(S3— 3 satellite range) ; and Region 3 near the 
equatorial plane (geosynchronous satellite meas- 
urements) . 



Region i 


Inverted V Precipitation 


The signature of an auroral arc in electron 
precipitation data is an inverted V* structure in 
an electron spectrogram ( Hoikkiln . 1970; Frank 
and Acker so n, 1971), In an inverted V the peak in 
the electron spectrum (monoenergetic peak) in- 
creases in energy toward the. arc center resulting 
in an Inverted V structure in an energy versus 
time spectrogram. Arcs or inverted V's are typi- 
cally a few hundred km in latitudinal width when 
detected aboard satellites, whereas rocket obser- 
vations, capable of better spatial resolution, 
reveal inverted V structures down to several km 
in width. There is some confusion in terminology 
in the literature concerning the identification 
of arcs with inverted V structures. Satellite 
papers refer to arcs imbedded in inverted V's, 
whereas from the perspective of rocket measure- 
ments any auroral form steady enough in time to 
traverse with a rocket payload has An inverted V 
structure. The difference appears to be one of 
measurement resolution. In fact, rockets have 
yet to identify the small scale structure (tens 
of meters) associated with optical arcs ( Maggs 
an d Da vis, £968), Finally, the relationship of 
electrostatic shocks seen by the S3-3 satellite 
( Hocor ct. al. , 1977) to inverted Vs or auroral 
arcs also remains unclear at the present time. 

Inverted V's occur on open as well as closed 
field lines. Figure 1 gives the latitudinal and 
longitudinal distribution of inverted V's seen by 
the AI3-D satellite ( Lin and Hoffman . 1979). They 
primarily occur in the evening sector auroral oval 
but are also present in the polar cap where field 
lines are. normally considered to be open, Using 
the definition that field lines are closed if 
they can support trapped radiation, inverted V's 
detected at satellite altitudes occur on both 
open and closed field lines ( Fran k and Ackerso n, 
1973; Mi a era et. al. , 1976; Eliasson ct. al. , 

1978; Lin and Ho ffma n, 1979). At rocket altitudes, 
however, inverted V's are generally seen coinci- 
dent with a mirroring population of energetic 
ions and/or electrons suggesting cloned field 
lines ( Venkatarngan et. , 1975; Whale n and 
Ofiley, 1979).' 

The monoenergetic peak in the. inverted V elec- 
tron spectrum is generally taken as evidence for 
electron acceleration by a quasi-static electric 
potential. There are. a number of observations, 
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however, that surest the inverted V electrons 
have passed through a turbulent region. The pitch 
angle, distribution of the monoenorgetie peak elec- 
trons is isotropic minus the loss cone, at low al- 
titudes ( Arnoldy ct. al . , 1974; Lin and Hoffman , 
1979). The temperature of precipitating electrons 
with energies above the spectral peak is propor- 
tional to the peak energy suggesting thi ; they 
are heated in the acceleration region ( Lin and 
Hoffman , 1979).. The monoenorgetie peak of in- 
verted V's represents a stable electron distribu- 
tion when observed at rocket altitudes ( Kaufmann 
et . al. , 1978) as the example in Figure 2 shows. 
Finally, the V-shaped pitch angle distributions 
observed by the AE-D satellite for energies below 
the monoenorgetie peak has been explained as due 
to a field-aligned component combined with iso- 
tropic backscattered electrons and a significant 
number of electrons with energies slightly below 
the peak energy trapped between their mirror 
uoints and an electric potential structure ( Lin 
an d H offman , 1979). 

The field-aligned electrons in an inverted V 
occur at energies below the monoenorgetie peak 
producing double peaked spectra as seen in Figure 
3 (Ar noldy e t . al, , 1974; L,in and Hoffman , 1979; 
Kaufm ann , 1980) . The 5 keV peak in the spectra 
of Figure 3_JLs isotropic In pitch angle, whereas 
the 2 •>• 3 keV peak is field -aligned and observed 
only at small pitch angles. This field-aligned 
component will be discussed in more detail later 
in this review. 

Ion Precipitation 

The. ESRO T data showed that ions and electrons 
simultaneously precipitated on the same field 
line, but l<eV ion flexes and spectra wore- less 
variable than those, of the electrons in the same 
energy range (Uultgvist , 1979). The ESRO t data 
also showed that keV electrons were frequently 
field-aligned in the. same region ns ions, however, 
more sporadic and narrower in space/ time than the 
ion field-aligned fluxes. Field-aligned keV ions 
occur at tliG edge of inverted V's with spectral 
softening at the center and a tendency toward 
isotropy (Hultqvtst, 1971). In contradiction to 
the ESRO I results cited above, the ESRO 4 and 
IH.ILF1 '5 data show no field-aligned ions and rarely 
any hot ions in inverted V's, respectively ( Frank 
and Acker son , 1972; ". Hasson et. al. , 1979). 

(1971) has interpreted the different 
results in terms of a solar cyclu effect. ESRO I, 
operating during solar maximum, observed both ions 
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electrons accelerated by a potential above the 
satellite. The potential accelerating the ions 
occurred below the election potential and was 
created by the precipitating clcctronn interact- 
ing with an elevated neutral atmosphere during 
solar maximum conditions. 

Rocket measurements of ion precipitation are 
more limited because of the low altitude obtained 
by many of the flights. Spectral and intensity 
changes are generally anticorrelated when ions 
and electrons precipitate simultaneously (Br yant 
et. al. , 1977; Edwards et. al. , 1978). A notable 
exception is the case of general ion ani electron 
anticorrelation until the center of an electron 
inverted V is reached where suddenly a band of 6 
keV ions appeared simultaneously with the elec- 
trons ( Bryant et. al, , 1977) . Ion data from a 
number of Canadian rocket flights ( Miller and 
Wha len , 1976) show relative independence of energy 
spectra from auroral activity. Pitch angle dis- 
tributions were observed to be either isotropic 
over the upper hemisphere or trapped, and there 
was, in general, little correlation with electron 
precipitation. 

Field-Aligned Precipitation 

As discussed above, the monoenergetic peak in 
the electron spectrum is generally isotropic in 
pitch angle. The electrons discussed here are 
sporadic bursts of intense field-aligned precipi- 
tation. The field-aligned electrons are usually 
less than a few keV in energy and are below the 
monoenorgetl.- peak in the spectrum ( Arnoldy et, 
al_,_, 1974; Lin and Hoffman , 1979). Figure 4 shows 
a flux ratio spectrogram of a satellite inverted 
V event. The heavy line tracing the spectral peak 
lies in a region of general pitch angle isotropy 
whereas the field-aligned fluxes are transient 
and of lower energy. The field-aligned fluxes 
generally occur over a range of energies (Arnoldy 
et. al . , 1974; Raitt and Sojka, 1977; Lin and 
Hoffma n. 1979). 

Field -aligned electrons arc generally seen dur- 
ing active, auroras (Mc Ewen . 1972; Arnoldy et. al.. , 
1974; Bosqued et. al. , 1974) and are localised on 
edges of arcs in regions between different plasmas 
(Wh alen and McDiarmid , 1972; Bryant et. al. , 1973; 
H aohlum and Moestue , 1973; Arnoldy , 1974; Lund in , 
1974; Pul Liam et. al ., 1980). Figure 5 gives two 
energy channels of a rocket electron detector 
along with ground photometer measurements in the 
lower panel. The rapid oscillations in the 1.4 
keV detector at 195 and 310 seconds are a spin 


modulation an the detector swept through the 
field-aligned precipitation localised near the 
northern and southern borders of the aurora. 

The field-aligned precipitation occurs in 
bursts or fluctuates at a few Herts ( McEwc n , 1971; 
A rno Id v e.t. al . , 1974; Bosqued et. nl. , 1974; Boyd 
and Davis, 1977; Rnitt and So.jka , 1977). If fluc- 
tuating, nil energies field-aligned fluctuate in 
phase (Li n and Hoffman , 1979). This later obser- 
vation rules out an energy dependent process for 
producing the field-alignment. Finally, the field- 
alignment is generally within 20° of the field 
line at rocket altitudes and on occasion less 
than 10° ( Bryant et . al . , 1973; Arnoldy et. al. , 
1974; Bosqued et. al. . 1974). To maintain such a 
degree of field-alignment the mechanism could not 
be too far removed from the measurements. This 
combined with their short-lived nature and field- 
alignment over a range of energies suggests that 
the phenomena is the result of the acceleration 
of ambient electrons in a large scale transient 
electric field. The relationship of the field- 
aligned fluxes to the inverted V acceleration re- 
gion remains unclear at the present time. 

Low Energy Electron Precipitation 

The electron differential directional intensi- 
ty belcw the monoenergetic peak in the precipita- 
ting electron spectrum can be fit to a power law 
as E“T where y varies between 1 and 2 ( Reason er 
and Chappell , 1973; Arnoldy and Cho y, 1973). Over 
an inverted V or a discrete auroral arc these low 
energy electrons arc generally isotropic over the 
full 180° range of pitch angles. Evans (1974) 
has interpreted this result in terms of an elec- 
tric field trapping secondaries and degraded pri- 
maries between the potential region and the at- 
mosphere. Recent efforts by Pulliam et. al . , 
(1980) to fit this model to rocket; data have been 
quite successful. 

Recent measurements over a diffuse aurora have 
revealed similar isotropic low energy electrons 
( Evans a n d Moore , 1979). Since the diffuse aurora 
electron precipitation has no peak in the spectrum 
tiie electrons apparently have not fallen through 
a pote.ntial region. There is, therefore, no mech- 
anism t.o reflect the. backscattered electrons at 
1 -*■ 2 earth radii above the aurora. To achieve 
the measured isotropic distributions, Eva ns and 
Moore (1979) have postulated that the backscat- 
terod electrons fill the entire flux tube and as 
a result of pitch angle scattering could be a 
source for the plasma sheet electron population. 



Ev ans and M oore (1979) have modeled the plasma 
sheet spcctrum*~and presumably the spectrum pre- 
cipitated during diffuse aurora are due. to a pri- 
mordial Maxwellian (Figure 6) and the atmospheric 
luckscattercd electrons that are in equilibrium 
with this Maxwellian. In the model spectrum the 
Maxwellian electrons dominate at hijgh energies 
and the atmospheric electrons dominate at low 
energies giving a spectrum which approximates that 
measured quite well. This topic will be discussed 
further under Region 3. 

When net over visible aurora the low energy 
precipitated electrons often remain but are not 
isotropic in pitch angle. Figure 7 gives down and 
up-going electron energy spectra measured by a 
sounding rocket in the region between two visible 
auroras. There is a definite downward streaming 
population of suprathermal electrons between 100 
and 500 eV energy which are not secondarily rela- 
ted to any higher energy component. Figure 8 gives 
data from another rocket flight where a very nar- 
row arc (spatially coincident with the 3 keV 
fluxes) was imbedded in widespread low energy pre- 
cipitation. Above the arc, at 145 seconds flight 
time, the electrons below 1 keV were isotropic in 
pitch angle while outside^the arc, at 200 seconds 
there was a net downward streaming of the less 
than 1 keV electrons. This net precipitation of 
less than 1 keV electrons carries significant 
current out of the atmosphere and could represent 
the high energy tail of the precipitation respon- 
sible for the large scale evening sector Region 1 
currents except when over an arc where the Region 
1 current is carried predominantly by the accel- 
erated energetic electrons. 

There is a growing body of evidence that supra- 
thermal electrons are rather easily heated by wave 
interactions. Measurements from both the Polar V 
and Echo III electron accelerator payloads give 
an intense background of suprathermal electrons 
that remained after the energetic beam injections. 
These electrons could not be interpreted in terms 
of beam secondaries nor as scattered and degraded 
beam primary electrons ( Maehlum et. a l., 1979; 
Arnold y and W in ckle r, 1980). The accelerator pay- 
loads are known, however, to be intense sources 
of wave turbulence when operated in the upper 
ionosphere (Ca rtwright and Kellog g, 1974). Such 
turbulence could heat the particle population near 
the payload. Johnso n and Sojka (1980) cite a case 
of ambient thermal electrons being heated to tem- 
peratures as high as 10 & K by the passage of a 
field -aligned beam of a few hundred eV through 
the ionosphere at altitudes over 660 km. A beam/ 



plasma interact: lan is cited as the. heating mech- 
anism. Finally, recent results we have obtained 
from the release of an argon plasma at 200 kin al- 
titude show considerable modification of the nat- 
urally occurring precipitation below 1 keV. Figure 
9 gives some preliminary electron measurements 
from this experiment. An isotropic population of 
electrons at a few hundred eV is introduced as 
well as a component below 100 eV streaming down 
the field lines. These electrons were closely 
correlated with plasma injections when wave in- 
struments aboard the payload also measured greatly 
enhanced intensities. Such modifications of the 
low energy component can alter the light emission 
and electrodynamic properties of the upper iono- 
pherc from that expected on the basis of atmos- 
pheric collisions above. 

Region 2 

The S3-3 spacecraft surveying the region of 
space from 400Q ■+ 8000 lem above the auroral arc 
has provided considerable, evidence for large 
scale electric, fields along the magnetic field 
from pitch angle and energy distributions of elec- 
trons and ions. The satellite is generally found 
to be immersed in an electric field ^ 1 volt/km 
extending oyj?,r a few thousand kilometers ( Chu 1 
and Sc hnl a, 1978; Clad is and Shar p, 1979). We now 
briefly review the different types of particle 
population measured by the S3-3 mission, not all 
of which are. readily interpretable in terms of 
quasi-steady field -aligned electric field. 

UF t (U p ward Flowing Ions ) 

UFI with a directional intensity J > 2 x 10 r> 
inn* 1 sec stor ko.V)“* are a common feature seen on 
approximately 60% of the S3-3 orbits. They coin- 
cide with the. statistical auroral oval and have a 
strong dawn-dusk asymmetry with a maximum in the. 
dusk sector (Ch ielmett i e t. al. . 1973). The pre- 
dominant species in the UFt are. 11* and ions 
(Sh elle y et. nl., 1976). UFI and electrostatic 
hydrogen cycle Iron waves coincide in over 90% of 
the events studied ( Kintnor et . al .. 1979). Fin- 
ally, UFI are observed to have a beam pitch angle 
distribution of (0 ■v 180°) or a conic pitch angle 
distribution (90° < 0 < 180°). 

UFI Beams 


UFI beams have a maximum frequency of occur- 
rence in the dusk sector and are observed prima- 



rily above .5000 km altitude ( Co rn ov e t.nl. , (1930). 
Beams are observed primarily in the region of up- 
ward f ie ld-ai igned current and coincide with elec- 
trostatic shocks in the evening sector (Cfittell 
ct. al .. 1979). It would appear that the beams 
are a phenomena associated with discrete auroral 
forms. 

UFI Conics 


Conics are observed at lower altitudes than 
beams with a tendency to prefer the higher alti- 
tudes during disturbed times ( Gorney et. al. , 

1980) . During magnetically quiet times conics 
tend to be a daytime cusp phenomena becoming more 
uniformly spread in local time as the activity 
increases ( Ghielmetti et. al ., 1979). Conics are 
observed in all field-aligned current regions 
( Cat tell et. al . , 1979) hence appear to be less 
associated with inverted V's than beams. Conics 
are apparently the same ions as the transversely 
accelerated ions (TAI) seen at lower altitudes 
( Whalen et. al. , 1978; Klumpar , 1979; Ungstrup 
et . al. , 1979), 

DFI (Downward Flowing Ions) 

* 

DFI occur much less frequently than do UFI in 
the aurora.l~"region and do rot display a dawn-dusk 
asymmetry ( Ghi e lmetti et, al ., 1979). DFI occur 
equatorward of UFI and in regions of Inward (eve- 
ning) and outward (morning) field-aligned currents 
( Cat tell et . al . , 1979), lienee appear to be asso- 
ciated with the diffuse aurora. DFI are correlat- 
ed with trapped keV ions and occur preferentially 
in spatially localized regions of enhanced hot 
keV plasma ions. Because of this, Ghielmetti, et 
al . (1979) hove associated DFI with the injection 
of ring current ions in the equatorial plane. Some 
DFI are associated with retarded electron preci- 
pitation suggesting an electric field accelera- 
tion opposite to that in inverted V's ( Ghielmetti 
et. q] . 1979). The first example of an ion in- 

verted V was presented at this conference by Bur- 
gess and Re if f (19S0) . 

Counte rs!: reaming Electron Beam s 

Counter-streaming electron beams have energies 
in the keV range. They are often spatially asso- 
ciated with regions of UFI conics ( Fennell ct. al, 
1979). They are field-aligned to within several 
degrees and are generally comparable in width for 
bot h upward and downward beams (Sha rp et .al 1980). 


Tim energy spectra for count au* stream i n 5 * electron 
beams are broad mid generally unstructured. In 
the region of eounterstruainlng electron beams, the 
electron beams, fc ha lows cone, of background elec- 
trons Is not affected, ruling out a quasi-static 
parallel electric field. Sharp ot. al. (1930) have 
suggested that the counterstreaming electron beams 
are locally accelerated ambient electrons by 
flickering double, layers. Counterstreamir.g elec- 
tron beams appear to be similar in many respects 
to the field-aligned electron bursts seen by low 
altitude satellites and rockets discussed above. 

Region 3 

Plasma injections at synchronous orbit are the 
equatorial plane signature of an auroral substorm. 
Although much correlative work between ground and 
synchronous satellite measurements has been done 
(A kasofu , 1977) , how one region maps to the other 
is still unclear. Moreover, the physics of the 
injection of particles at synchronous orbit re- 
mains to be. defined. That there is an exchange 
of particles between the equatorial plane and the 
ionosphere is certain, however, the degree to 
which the ionosphere provides plasma to the equa- 
torial plane (plasma sheet) is yet to be fully 
explored. As we have seen, measurements at 'v 1 
Rr (Region 2) reveal many particles leaving the 
ionosphere. Presumably these particles are those 
that have been detected at synchronous altitude 
as "source cone" ions and electrons ( Me 111 wain , 
1974 ; M auk and Mclllwnin , 1975). It is to be 
noted, however, that the source cone, particles 
observed have a wider pitch angle spread than the 
actual loss cone at 6.6 Rg. If the particles 
were accelerated in Region 2 they must have been 
diffused in pitch angle, before reaching the 
equator. 

An example of an electron injection at syn- 
chronous orbit is given in Figure 10. At 04 12 UT 
the intensity of the higher energy electrons in- 
creased by one or two orders of magnitude while 
the low energy electrons (< few kcV) decreased by 
about as much. The injection represented a drama- 
tic change in spectrum from a monotonically de- 
creasing one to nearly a Maxwellian. Following 
the injection, die spectrum then gradually *e- 
turned to approximately its old shape and to pre- 
injection intensity levels. If an aurora wag 
initiated as a result of this injection, then the 
gradual return of low energy electrons might re- 
present the filling of the plasma sheet with back- 
scattered and secondary electrons from the auroral 
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ionosphere ns suggested by 1 lv; n w > ai u 1 ! 1< x ) r e (1979). 
Concluding Remarks 

This review of auroral particle precipitation 
raises three questions about auroral arc forma- 
tion, First, what is the role of AC phenomena in 
auroral particle acceleration? Examples of bursts 
of field -aligned electrons, simultaneous ion and 
electron precipitation, counter streaming elec- 
trons, and conics (transversely accelerated ions) 
all suggest AC phenomena. Second, what is the 
importance of electrons < 1 kcV in the auroral 
beam? Clearly they are not all secondaries. In- 
verted V's appear to be imbedded in and possibly 
develop as a result of this precipitation. Final- 
ly, what role does the lower ionosphere have on 
the. formation of arcs? Recent controlled experi- 
ments with sounding rockets which have altered 
the lower ionospheric conductivity often have re- 
sulted in some form of anomalous particle preci- 
pitation. In addition, there is a growing evi- 
dence for fairly low altitude acceleration of 
particles occurring naturally, such as bursts of 
low energy field-aligned electrons and transverse- 
ly accelerated ions. 
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Figure Captions 

Fi e,. 1. Spatial occurrence map of 280 inverted 
V events. Ixscal time, and circles of constant 
invariant latitude are shown as reference (from 
Lin and Hoffman, 1979). 

Ili* Example to show how a peak in the dif- 
ferencial energy spectrum is not necessarily an 
unstable structure in the distribution function 
f . 

Sequential electron spectra from 
flight 18:91. Pitch angle of the electrons is 
given by the smooth curve. The lower two dec- 
ades of the flux scale are used to give pitch 
angles from 1° to 100° (from Arnoldy et. al. , 
1974) . 

Fig . 4. Spectrogram displaying the flux ratios 

J(-7°)/J(60 < ') of an inverted V event. The flux 
ratios are separated into four categories as 
shown in the upper right corner (from Lin and 
Hoffman , 1979). 

Fuj_. 5_._ Top two panels, directional intensi- 

ties 1.4 and 4.5 koV electrons. Bottom panel, 
ground meridian-scanning photometer looking at 
font of rocket field lino, energy of peak cur- 
rent (from Arnoldy, 1977). 

Measured electron spectrum during a 
diffuse aurora together with a model spectrum 



which is the sun of a JlnsvcU i.i« ami the dia- 
tribe! ion of atmospherically produced back- 
seat t need and secondary electrons that in in 
equilibrium with thin Maxwellian. The energy 
spectrum for the Maxwellian of temperature 2.5 
keV and density 0,2 cm”* is also plotted for 
comparison (from Evans and Moore, *1979). 

Spectra of electrons moving up and 
dovarThe magnetic field line measured simulta- 
neously by rocket detectors. 

Fig. 8, Rocket measurements of the precipita- 
tion Vf 0.3 and 3 ke.V electrons and simultane- 
ous up and down spectra taken at two times 
during the flight. Visible aurora coincided 
only with the 3 keV precipitation bursts cen- 
tered on 150 seconds flight time. 

£ig. r _J)_. Background auroral spectrum (heavy 
Tine) "and three simultaneous spectra taken by 
detectors scanning the pitch angle, ranges In- 
dicated during an argon plasma release. The 
auroral spectrum is dramatically altered below 
] keV. 

Fig. 10. ATS-d electron spectra as a function 
of time during a substorm injection event. The 
top eight traces give count rates for eight 
energy channels as a function of time. The 
energy in keV is given to the right of each 
trace. — 
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Modulation of Terrestrial Ion Escape Flux Composition (by Low- Altitude 
Acceleration and Charge Exchange Chemistry) 

Thomas E. Moore 

Space Science Center, University of New Hampshire, Durham, New Hampshire 03824 

Motivated by recent observations of highly variable hot plasma composition in the magnetosphere, 
control of the ionospheric escape (lux composition by low-altitude particle dynamics and ion chemistry 
has been investigated for an e , If, 0 + ionosphere. Theoretical polar wind results of Banks and Holzer 
and observations of highly nonlhermal ion escape have been used as guides. It is found that the fraction 
of the steady state escape flux which is O* can be controlled very sensitively by the occurrence of parallel 
or transverse ion acceleration at altitudes below the altitude where the neutral oxygen density falls rap- 
idly below the neutral hydrogen density and the ionospheric source of O' tends to be rapidly converted 
by charge exchange to H'. The acceleration is required both to overcome the gravitational confinement 
of O' and to violate charge exchange equilibrium so that the neutral hydrogen atmosphere appears ‘opti- 
cally' thin to escaping O'. Constraints are placed on the acceleration processes, and it is shown that O' 
escape is facilitated by observed ionospheric responses to magnetic activity. 


INTRODUCTION 

The hot magnetospheric plasma has a variable composition 
which usually includes and is at times dominated by heavy 
terrestrial ions, most probably O', Yet models of ionospheric 
escape (e.g., Banks and Holzer [1969b], hereafter referred to as 
BH) predict that only a fraction of order 10 J or less of this es- 
caping flux should be O'. This is in accord with a comparison 
of ionospheric thermal speeds and the gravitational escape ve- 
locity. Though auroral zone particle acceleration delivers 
energies much larger than either thermal or gravitational 
binding energies, BH have shown that there is another ob- 
stacle to O' escape, namely, the charge exchange O' flowing 
up from its source on the neutral H atmosphere, This paper 
presents an analysis in which ion acceleration is regarded as 
an imposed free parameter. Steady stale conditions are identi- 
fied under which the ionosphere will respond by releasing pri- 
marily O' rather than IT . 

The compositional variability of the hot magnetospheric 
plasma is illustrated by observations from geosynchronous or- 
bit [Ceto et al„ 1978], polar low-altitude orbit [Shelley et al„ 
1972; Gliielmetti et al„ 1978], and rocket altitudes [Lynch et al,, 
1977; Moore and Evans, 1979]. These observations include sit- 
uations with O'/H' ranging from less than a few percent to 
greater than unity. The correlation of O' content with mag- 
netic activity reported by Sharp et al. [1976] suggests than an 
ionospheric response to magnetic storms may be involved. 
Moreover, unusual enrichments ofN' and molecular ions ob- 
served in the thermal ionosphere [ Hoffman et al., 1974; Taylor 
et al., 1975] during magnetic storms further support such a hy- 
pothesis. 

The conceptual separation of energization mechanism and 
compositional control is supported by observations of keV 
ions jetting up out of the ionosphere between 1000- and 8000- 
km altitude above the auroral zone [Kintneret al., 1979; Klum- 
par, 1979; Sharp et al., 1977; Shelley et al., 1976] from a vaiiety 
of polar orbiicrs. These ions are generally field aligned to a 
high degree, sometimes showing ‘conical* pitch angle distribu- 
tions with minima at 90° and 180°, The point of most impor- 
tance here is that the ion jet composition is highly variable 
from essentially pure H' to pure O' (see especially Kintner et 
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al. [1979]), Unp strap et al. [1979] have proposed that these ion 
flows result from transverse acceleration of ions below the sat- 
ellite by ion cyclotron waves driven unstable by field-aligned 
currents. The resulting magnetic mirror force then ejects the 
ion upward out of the ionosphere. While gyro-resonant pref- 
erential acceleration of one or another ion species is an attrac- 
tive model for the observed composition variations, such ef- 
fects are ignored in the present treatment, which presumes 
that acceleration is effective for al! species. Instead, attention 
is focused on ion charge exchange chemistry, with results 
which have implications for gyro-resonant effects as well. 

The approach has been inspired by the results of BH, which 
represent a solution to the problem of ion dynamics and 
chemistry when no acceleration is imposed beyond the natu- 
ral thermal escape of ionospheric plasma. Those results are re- 
viewed in the next section. In subsequent sections, continuity 
of O' is considered for the more general case of an arbitrarily 
imposed vertical flow profile, and the results are evaluated for 
special cases of (1) parallel acceleration and (2) transverse ac- 
celeration. A final section summarizes conclusions which may 
be drawn from the results obtained. 

Polar Wind Theory (High-Latitude 
Vertical plasma Transport) 

Considerable theoretical work has been done on the escape 
of ionospheric plasma into the high-latitude magnetosphere 
beginning with, for example, Banks and Holzer [1969a], and 
the escape phenomenon has been termed the polar wind by 
Axford [1968] in recognition of a certain analogy with the so- 
lar wind. Lemaire and Scherer [1973] have reviewed the lively 
debate which followed the initial papers. This debate centered 
on the relative merits of hydrodynamic and kinetic ap- 
proaches to the problem, and also exhibited some analogies 
with catlief debate surrounding solar wind theory. Hydro- 
dynamic models are evidently appropriate at low altitudes 
where collisions are important, while a kinetic approach is ap- 
propriate in collisionless regions, the two approaches being of 
course complementary rather than contradictory. Certain 
quantitative refinements involving ionization rates and neu- 
tral atmosphere models have recently been introduced by 
Raitt et al. [1978], with particular attention to helium escape, 
The emphasis of the present paper will be on H' and O' es- 
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cape and the accidentally resonant charge exchange chemistry 
which couples these two ions and their neutral atoms. 

The hydrodynamic approach has been applied to a polar 
ionosphere of H + , Hc\ O', and e by BH, who performed a 
steady state calculation using continuity and momentum 
equations for each species. The equations used describe verti- 
cal motions taking account of sources and losses of ionization, 
charge exchange chemistry, and ion-ion, ion-neutral friction. 
Electron inertia is ignored, and ions and electrons are as- 
sumed isothermal except that separate electron and ion tem- 
peratures are specified. Pressure is taken to be isotropic and to 
obey a simple gas law. The neutral atmosphere is described 
with the model of M. Nicolet ct al. (unpublished manuscript, 
1969). A dipolar B dependence is used to describe the varia- 
tion of flux tube cross section with altitude, and the plasma 
pressure at great altitudes is assumed to be sufficiently low as 
to force n transonic solution for the light species H + , Hc + , The 
result of the calculation is a specification of and Uj(z\ for 
z up to ~ 1 and where j enumerates all species. 

The results of their calculation can be summarized as fol- 
lows: A large escape flux of H + occurs amounting to over 
10*H '/cm 2 s. Escape fluxes of Hc + and O* are much smaller 
(by 10' a and 10 \ respectively) Most importantly for the 
present work, a large upward flux of O* (>10*/cm* s) occurs 
at altitudes below ~6Q0 km which, through charge exchange 
on neutral H, supplies the escaping H + flux. It may reasonably 
be inferred that the ionosphere has a built-in obstacle to the 
escape of O* in the form of a neutral hydrogen component 
which tends to neutralize any upward flowing 0 \ This result 
tends to make the observations of escaping O* cited earlier 
very difficult to understand. 

Wherever the vertical flows are slow enough (much less 
than the smallest scale height per charge exchange period), 
charge exchange equilibrium will be maintained, implying 
that [BH] 

n(H*)„ 9n(H)(T n V'* 

7(0 + ) " 8«(0) \ T,j 

where the n are number densities and T„, T t are neutral and 
ion temperatures. This condition is satisfied in the BH model 
below ~60G km, where n(0) drops below «(H) (the ‘H, O 
crossover level’) and 0 + is converted to IT, One might sus- 
pect that a lower 11 density might shift the equilibrium in fa- 
vor of 0 + , but this is not the case. It must be borne in mind 
that the polar wind H + escape is driven by the thermal energy 
of the ions (though it is crucially aided by the 0 + -e" gravita- 
tional separation E field). Normal thermal escape of 0 + is 
very small. BH have modeled the low H case by considering a 
higher than normal temperature neutral atmosphere. Figure- 1 
shows model neutral H atmospheres used by BH and best cur- 
rent models from Tinsley [1974], with neutral O crossover lev- 
els indicated by dots. A modest temperature increase can 
clearly produce a dramatic decrease in H above 250 km, due 
to increased neutral H escape. The result is a reduced source 
and hence escape of H + . This reduces the demand for and 
hence the flux of 0 + , The charge exchange region moves up- 
ward to the new n(H), and «(0) crossover level, and 0 + es- 
cape is still small, The combination of gravitational and 
chemical obstacles would appear to effectively prevent steady 
0 + escape, regardless of forces imposed at altitudes above the 
n(0) ~n(H) level where IT is produced. 

However, it seems reasonable to inquire whether the iono- 



Fig. I. Model neutral hydrogen atmospheres for altitudes above 
200 km, for three temperatures. Profiles from Banks and llolzer 
[196%) are marked BH, Those from Tinsley [1974] are marked T, 
Rough scale heights are indicated on profiles which are nearly straight 
lines, and the altitudes at which neutral 0 (scale height ~70 km) falls 
below neutral H are indicated with dots. 

spheric escape can be driven out of charge exchange equilib- 
rium so as to favor 0 + escape by application of acceleration at 
altitudes where the 0 + upflux exists. Auroral acceleration is 
known to impart energies far in excess of ionospheric thermal 
energies, thus violating many of the conditions for hydro- 
dynamic flow. As a means of investigau tg low-altitude accel- 
eration effects, the equation of motion will be ignored and the 
dynamics regarded as a. free parameter, Only a simplified con- 
tinuity equation for 0+ will be considered. 

Ion Escape From a charge Exchanging 
Atmosphere 

In this section the ion dynamics will be parameterized in 
terms of a specified vertical velocity profile v,(z). A continuity 
equation for some initial upward flux of 0 + under the action 
Of forces appropriate to produce v,(z) will be written including 
the effects of charge exchange loss and gain due to collisions 
of the 0 + with H and the fast O thus produced with ambient 
H + , Later, the solution obtained to this equation will be eval- 
uated for some simple assumed upward velocity profiles. 

The steady state continuity equation for vertical ion trans- 
port in a plane parallel atmosphere is appropriate here in view 
of the long times over which 0 + is observed to vary in the 
magnetosphere and the fact that most of the charge exchange 
occurs below 1000 km. It is 



where J 0 * « rio-v,. It is assumed that the z axis is parallel to 
the magnetic field, which will be assumed essentially vertical, 
thus restricting this analysis to a high-latitude situation. A 
term whicn accounts for field line divergence with altitude has 
been ignored here, since its effect is to account for geometri- 
cally produced changes in J Q + which are not of concern here. 
All sources and sinks of 0 + other than charge exchange will 
be taken to be in balance everywhere. Neutral hydrogen is 
present with number density protons with density n„+, and 
the cross sections for the forward and reverse reactions 0 + + 
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H » 6 + H* are taken as a and o„ respectively. Continuity 
then demands 

^7*- " •“•Hh<<7V}«o* + M,|»<0,V>rt 0 


more 0 * survives, the main source of H* is removed, and the 
Q* source term due to the reverse reaction becomes even 
smaller. For small w the solution for 0’ survival when gravi- 
tational escape is assured by external forces is just the fraction 
of initial 0 * flux remaining after charge exchange loss: 


Only the fast 0 produced by charge eachange of the initial 0 + 
flux (/, « («„,»’,),) can become & source of O* flux by reverse 
charge exchange. Hence the appropriate oxygen number den- 
sity for the third (source) term is given by 

(Wo v t) ““ Jf ~ (Wj)*Vi) 

The factor (ov) indicates the product of a nnd the root square 
sum of ion thermal, neutral thermal, and ion-neutral trans- 
lation speeds. 

Defining /+ m J 0 ,/J„ the equation for continuity becomes 

v> iL m ~„ M <ov>/ + + n„.(o,v>(l - /♦) 

Solution of this radiative transfer type equation is simplified if 
one defines a charge exchange depth by analogy with optical 
depth: 


T * 




dz 


In radiative transfer the ratio v ,/v is usually taken out of the 
integral as a cosine factor. However, since it is desired to allow 
for an ion pitch angle which may vary with z, this factor will 
be retained in the integral here. Note that r represents an at- 
tenuation depth measured along the path of an O* ion in es- 
caping from an initial height z,, to large altitudes. With an ad- 
ditional definition, 


— 2a 4=4 

n,,(ov> 

the continuity equation becomes 
d!+ 

~~~ = -/+(1 + w ) + H> 

(IT 

Here w is properly viewed as a function of altitude. A closed 
form solution does exist and may be integrable for certain 
simple forms of w(z). Proper calculation of w requires consid- 
eration of the continuity and dynamics of H + , which in the 
calculation of BH had its major source in the 0 + -H charge ex- 
change, It can be argued from their results that upward escape 
of H + occurs so rapidly that W densities remain very low and 
w never becomes large enough to influence the final result for 
0 + . This may be seen as follows: In the case w - constant, the 
solution for / + (Q) “ 1 is 

r+ — !L_ 1 1 + J_ e (H»)t 

1 1 + H-\ 1 + H- 

This in fact is not a bad approximation, since H + and H have 
quite similar altitude distributions above 400 km in the BH re- 
sults. The solution has the appropriate properties for small 
and large t, i.e„ / + ( 0) = 1 and / + (°°) — » h'/( 1 + w). Small f* 
at high altitudes, as found for the polar wind, carries the im- 
plication then of small w\ In fact, n (( « may be seen to be much 
smaller than «„ in the BH results. The source term in the con- 
tinuity equation is evidently kept small by rapid upward loss 
of H + , even when no external force assists in H + transport. As 


I* » e r 

Here, r is the charge exchange depth encountered between 
some starting altitude and a higher obseivation point, 
Clearly, the amount of O* escaping from the atmosphere un- 
der the described dynamical conditions is sensitively depen- 
dent on t. It remains to determine how r may vary in response 
to hypothetical dynamic situations. This requires assumptions 
about the form of v.(z). 


Parallel Acceleration 


The simplest possible assumption for v,(z) is that it is a con- 
stant which is imposed by some initial impulsive acceleration 
and maintained against gravity by some continuing force. 
This is a very good approximation to a situation including a 
low-altitude plasma double layer and a useful idealization of a 
situation involving a distributed ver'ical potential drop which 
extends to low altitude, In such a case, 


r 




ii 


with A r n the column density of neutral hydrogen above the ac- 
celeration lower edge. In terms of an assumed exponential at- 
mosphere /in ■* n 0 e *"■'»; t ** z 0 /z„ where z, * [«ov)/v,)n 0 ] 1 
defines a ‘charge exchange survival height,’ Appreciable 0 + 
escape then requires z, -a z,„ i,c„ 1, That such is not the 
case in the calculation of BH may be seen by substitution of 
appropriate values: N u ~ 10 M («» ~ 1,4 X 10 6 , z 0 ~ 750 km), 
v, ~ 2 X 10 4 cm/s, (av) - 4.3 X 10 " (7*„ *F 7V/I6 + 1,2 X 
10 ’V, 2 )" 2 - i ,3 x 10 v with T„ *=750°K, T 0 * - 3000°K. The 
result for z, is —110 km yielding r » 6,6 and an 0 + survival 
fraction ~ 1 Q \ in rough accord with the complete physical 
model. 

This result indicates the conditions under which an appre- 
ciable fraction of the terrestrial ion escape flux can be 0 + . 
First, an externally imposed accelerating force must overcome 
the gravitational confinement (10 eV is required for 0 + ), 
Next, the charge exchange depth of H must be reduced to a 
state of ‘optical thinness,* This can be accomplished in one or 
both of two ways: ( 1 ) a reduction of the neutral hydrogen col- 
umn density above the H, 0 crossover altitude or (2) a reduc- 
tion of (av)/v, by acceleration at or below that level. 

The first of these possibilities involves an alteration of the 
neutral hydrogen atmosphere. Tinsley [1974] has reviewed 
current knowledge of the hydrogen ‘geocorona’ and its varia- 
bility, Polar wind H + loss may in time deplete neutral H, The- 
oretical models predict that the column density above 250 km 
is tremendously variable with (isothermal) temperature due to 
escape of fast neutrals (refer to Figure 1). Drops of an order of 
magnitude result from a temperature rise by only a factor of 2 , 
Storm activity is expected to dominate the auroral ionosphere 
energy budget producing strong heating at low altitudes. 
However, Tinsley notes optical measurements indicating max- 
imum column density reduction below 450 km of ~40%. 
Heating and escape of topside hydrogen may be washed out 
by flows from other altitudes and latitudes. In any case, with- 
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out imposed demand for ions (upward acceleration), a reduc- 
tion of N tl simply reduces the IP escape flux and hence de- 
mand for O', so that no increase in O' escape results. 

The second possibility is more amenable to analysis. It 
could be effected by a dependence of 0 on v or an increase in 
v 4 . Figure 2 has been constructed using the BH expression 
given above and cross-section measurements for keV O' cited 
by Tinsley [I976J as indicated. In Figure 2, (ov> has been 
plotted versus the energy with which O' approaches cotd neu- 
tral H target atoms, though in fact neutral and ion thermal 
speeds must be root square summed with ion translation 
speed. The BH expression assumes 0 is constant with energy, 
and the keV measurements suggest that a has at most a weak 
energy dependence (~p or is w0 "). An expression for t(v,) 
may now be written, which is for the BH cross section: 

w 4.3 X IQ' "(T h > 4- 7W16 + 1,2 X 10' V) 1 '* 


This and the corresponding relation for the keV measurement 
extrapolation have been plotted in Figure 3 versus the O' en- 
ergy corresponding to v,. In both cases, JV„ has been taken as 
10 M cm s . Both curves decrease with B, for E t below the ther- 
mal energy of -0.1 eV, where they reach a value near unity 
and begin to saturate between unity and 0 . 1 . Increased ther- 
mal energies have the effect of raising the left end of the 
curves while leaving the right ends stationary, as indicated, 
The entire curves move upward or downward in proportion to 
the hydrogen column density N, as indicated at center. The 
curves demonstrate that the ionosphere considered by BH will 
appear ’optically' thin to O' ions escaping with translation en- 
ergy comparable to or larger than their thermal energy, i.c„ 
those escaping with small pitch angles. Hence an acceleration 
mechanism which can produce supersonic O' outflow at alti- 
tudes below the crossover of neutral 0 and H densities can to- 
tally alter the composition of the escaping ion flux, A surpris- 
ing aspect of this result is that only a relatively small upward 
energization (by auroral standards) is needed, provided con- 
tinued forcing at least maintains the outflow against gravita- 
tion, 

The most obvious mechanism for such a process is the mag- 
netic field aligned (parallel) electric field thought responsible 
for downwai .1 electron acceleration [Evans, 1974). However, 
the distribution of such parallel electric fields (£„) at low alti- 
tudes is subject to great uncertainty due to a lack of either ob- 



Fig. 2. Charge exchange cross sections for the reaction of 0 + + IT 
— * O + H' as a function of energy with which O' approaches H. The 
upper line corresponds to an extrapolation of low-energy values un- 
der the assumption of constant o. The lower curve is an extrapolation 
of measurements given by Tinsley [1976] for keV energies, suggesting 
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Fig, 3. Charge exchange depth of the hydrogen atmosphere at a 
column density depth of 10 M cm 1 ns a function of O' energy for two 
types of ion acceleration. E is the upward translation energy in the 
parallel acceleration case, the initial transverse energy imparted to the 
ion in the transverse acceleration case. 


servations or models of their behavior below ~2000 km. 
Clearly, a very small fraction of the inferred several kV poten- 
tial drop can be of importance, It is in the altitude range of in- 
terest (200 2000 km) that vertical current conduction makes a 
transition from resistive to collisionless behavior. If the is 
associated with the collisionless regime, penetration of E t to 
low altitudes may be enhanced by depletions of topside ioni- 
zation and/or neutrals, Plasma depletions have been observed 
at auroral and subauroral latitudes [Taylor et a!., 1975; Rao ei 
ai, 1978; Rob to el al„ 1971], notably in association with mag- 
netic activity. Further observational and theoretical work on 
the variability and conduction properties of ihls transition re- 
gion is clearly needed to assess the importance of B 5 pene- 
tration down to charge exchange altitudes. 

Transverse Acceleration 

It is of great interest whether the transverse acceleration by 
ion cyclotron waves proposed to account for observed conical 
ion pitch angle distributions can also successfully overcome 
gravitational and charge exchange barriers hindering O' re- 
moval. It is immediately clear that (1) again, the acceleration 
must occur below the H, O crossover level to have access to 
the low-altitude O' source; ( 2 ) the acceleration must occur in 
a lime much shorter than the local charge exchange collision 
period l/i|i(o r v >]~ 1 ~ 50 s for BH parameters. Next, the ion 
charge exchange depth must be evaluated for O' ions escap- 
ing under the influence of the mngnctic mirror force and grav- 
itation. With the following assumptions, (1) cla/elE = 0, 
(2) «h itoe" m \ (3) B is parallel to g and z, (4) there is dipole 
|fl| variation with altitude, and (5) ions are rapidly accelerated 
to E xa and have vanishing initial vertical bulk flow, the result 
for r (see appendix for derivation) is 

T " aN "i\-UEjE x jr 

valid for /; » 1, where p & R F Jz a and is expected to be of or- 
der 6,5. E 10 is the initial transverse energy, E t is the O' escape 
energy ~1Q eV, and 0 , N H have their usual meaning. This 
function has a singularity for £ i0 *=> y £,, due to the impossibi- 
lity of gravitational escape for E xa at or below that level. The 
behavior is plotted versus E in Figure 3, For larger energies 
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(he charge exchange depth drops abruptly to approximately 
unity for BH parameters, where it levels out. This mechanism 
thus , sutlers a disadvantage by a factor of order 2.5 compared 
with the parallel acceleration case in the limit of large acceler- 
ate , Moreover, a larger threshold exists for transverse accel- 
eration, since the fttll escape energy must be provided at once. 
Heating of 0* at rates slower than 10 cV/min merely hastens 
their loss by charge exchange, producing fast neutral 0, The 
fast neutral O may transfer energy to neutral H and thus alter 
its distribution, a$ discussed above. However, with heating 
rates much larger than 10 eV/min, occurring at altitudes be- 
low the H, 0 crossover level, the magnetic mirror force can 
convert an appreciable fraction of the transverse energy to up- 
ward translation so that the charge exchange depth drops to 
approximately unity and 0* escapes without much loss. 

One problem with this mechanism is the low altitude at 
which it must operate for the ionospheric parameters of BH 
(H, O crossover at ~600 km). Kindel and Kennel [1971] have 
found that the critical current density for ion cyclotron wave 
growth reaches very high values of 10 '' to 10 2 A/m 1 at 400 
km. They identified ionospheric variability as crucial to the 
minimum altitude at which waves become .-destabilized by a 
given current density. The waves are unstab s when a critical 
clcctron-ion drift is exceeded, i,e„ above an altitude which de- 
creases with decreasing electron number density (or increasing 
current density). Their results indicate that wave growth at 
400 km can occur only for the extreme case of an ionospheric 
trough condition with a current of 160 /rA/m 1 . Such condi- 
tions may occur during major magnetic storms. However, it 
still remains to show whether ion transverse heating rates are 
sufficient to reduce the 0* charge exchange depth, a task 
which is beyond the scope of this paper, Heating rates of 10- 
1507yion/s obtained by Drummond and Rosenblutb [1962], 
Rynn el al, [1974], and Palmadesso el al. [1974] and cited by 
Ungstrup e, al, [1979] appear to be ample; and the observa- 
tions of IVhalen el al, [1978] indicate that such heating can oc- 
cur at 400-600 km in the expansive phase of a substorm. 

Conclusions 

The notion that high-latitude ionospheric escape flux com- 
position may be influenced by the interaction of ion dynamics 
and charge exchange chemistry has been explored for an H\ 
0\ e ionosphere. The polar wind model of Banks and Holier 
[19696] for the thermal escape of plasma and observations of 
very nonthermal ion escape have been used as guides. A 
simple description of ion continuity with arbitrarily imposed 
dynamics has been used to argue that 

1, Appreciable steady escape of 0 + from the ionosphere 
requires ion acceleration at altitudes below that where the 
neutral O density falls rapidly below the neutral H density 
and the ionospheric 0 + source is converted by charge ex- 
change to H\ This level is normally in the 500- to 1000-km 
range, according to atmospheric models. The acceleration is 
required both to overcome the gravitational confinement of 
O* and to violate charge exchange equilibrium so that the 
neutral H atmosphere appears ‘optically’ thin to escaping 0 + . 

2, Parallel vertical acceleration to sO.l eV will suffice if 
continued forcing insures gravitational escape. Such contin- 
ued forcing is automatic if a kilovolt parallel potential drop 
exists, 

3, Transverse acceleration must occur at a rate much 
larger than 0.2 eV/ion/s in order that gravitational and 
charge exchange barriers be overcome. Rates still larger than 


this will not achieve much additional increase in O f escape. 
The transverse acceleration must stop as soon as the ion be- 
gins to move rapidly upward. 

4. Depletion of the neutral H atmosphere by heating and/ 
or escape my play a crucial role by (a) raising the H, O cross- 
over level so that O* persists to higher altitudes where acceler- 
ation mechanisms are known to be operative and (b) decreas- 
ing the H column density above the level where acceleration 
begins, 

5. Depletion of the topside thermal plasma may alterna- 
tively serve as an 0 + ‘throttle’ by lowering the region where 
acceleration occurs. 

This analysis shows that ionospheric charge exchange 
chemistry is an important factor which can effect strong mod- 
ulation of ionospheric escape flux composition by interaction 
with low-altitude acceleration processes. Resonant preferen- 
tial acceleration of a particular ion species need not occur for 
this modulation to be effective, and this mechanism may mask 
effects of such selective acceleration by consuming the accel- 
erated ions to produce fast neutrals. It can easily be appreci- 
ated that an analysis of a more complex ionosphere contain- 
ing NO* or other ions should reveal other interesting effects 
due to changes of chemistry with varying dynamics. These re- 
sults point to a bona fide interaction between magnetosphere 
and ionosphere which can potentially influence the oxygen 
chemistry of the ionosphere over storm or solar cycle time 
scales and that of the atmosphere over geologic time scales. 
Determination of the relative importance of plasma and neu- 
tral depletion or of low-altitude acceleration requires a better 
understanding of ionospheric response to magnetic activity 
and of parallel current conduction in the topside ionosphere. 
Additional observations of dc fields, waves, and particles in 
the altitude range from 500 to 1500 km, i.c., above current 
rocket or radar altitudes and below most satellite altitudes, 
would be most helpful. A desirable theoretical goal would be 
an extension of polar wind theory to situations with vertical 
currents and particle acceleration to well above thermal 
energies. 


Appendix: Charge Exchange Depth por 
Transversely Accelerated Ions 

Beginning with the definition for charge exchange depth, 




and expressing in terms cos 0, where i is taken parallel to B so 
0 is the ion pitch angle; and n, the ‘target atom’ number den- 
sity, is taken to be n •* n„e we have, with the change of 
variable l = z/z„, 


f oo 

r •* aN / dl c '/cos 0 
Jo 


Now cos 0 = [1 + (Vj. 2 /r/)]~ l/2 . A dipolar field is appropriate 
at altitudes of interest, so that 



Here, r 0 is the geocentric radius at which ions are accelerated 
to initial transverse speed v l0 . To get v, 2 , we assume that B - g 
= 1 and write an effective potential for field-aligned motion 
under the influence of gravity and the magnetic mirror force. 
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Here, g 0 is the acceleration of gravity at the starting radius r 0 , 
and 0» is chosen such that (/>,«■ 0 there! 






At r > r, )t B, + <}>, m 0 expresses vertical energy conservation 
f,v ions with v, ** 0 at r u . Thus 


m 






and 


cos O' 


I + 


Via 8 (ro/rY 




-2gv<ill--(r„/r)] 


with the substitutions r *• r 0 + z, t « z/r„, l> “ t 0 /z 0 , so that z is 
measured from r«, some manipulation of the expression for t 
gives 


T 



(1 + i/p) 3 

l-(2W' , io i )(i“P/0 ' 


Expanding this for small t/p and retaining terms to first order 
only, one obtains 


ruaN 

valid for p large, or 
, P hi 

73U - W3£jo)] 




‘’xo 2 IP. 


r at aN * 0/V : 


P 


1/3 


[I = {2E f /3E X0 )Y 


where E, * 1/2 mv? » tng„r 0 is the gravitational escape en- 
ergy from radius r 0 i E 1(l is the initial transverse ion energy, 
and N is the column density of target atoms above r„; p is the 
ratio of r a to the target atom scale height and is expected to be 
larger than 5. 
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A Nikc-Tomaliawk rocket was launched north over two quiet, late evening auroral arcs in March 1975. 
A northward magnetic disturbance was observed on the ground under the rocket trajectory. South of the 
arcs the northward electric field was 60 mV/m, indicating strong westward plasma flow. An eastward 
electrojet current layer was penetrated in the upward (light. Precipitating electrons were observed over 
each arc. The electric field decreased to below 20 mV/m over each arc and recovered to 50 mV/nt be- 
tween the arcs. Using the observed electron flux and a model of the ionosphere, the Hall and Pedersen 
conductivities were calculated. These conductivities were used, with the observed electric field, to calcu- 
late the eastward and northward components of the horizontal ionospheric currents. The eastward cur- 
rent calculated south of the first arc agreed well with the observed eastward electrojct current. The power 
dissipated by the Pedersen current S p |E| J was also calculated and compared with the power carried by 
the precipitating electrons. The Joule power decreased abruptly in tire auroral arcs, as the precipitating 
electron power increased. The total dissipated power wns approximately the same inside the arcs, be- 
tween them, and south of the luminosity. North of the aurora the electric field and dissipated power re- 
mained low. Field-aligned currents carried by the observed electrons were about a factor of 3 lower than 
those inferred from the magnetic field measurements, Likewise, current continuity arguments to keep the 
auroral current system divergence free required a larger field-aligned current than that obtained from the 
particle detectors. 


INTRODUCTION 

Plasma in the ionosphere drifts toward the west in the eve- 
ning auroral oval under the influence of a northward electric 
field, In the lower ionosphere, near 1 10-kin altitude, where the 
positive ions are retarded by collisions, the westward electron 
motion constitutes an eastward ionospheric current, A region 
of eastward drifting plasma is sometimes observed, to the 
north of the westward drill, with an associated westward iono- 
spheric current [Maynard, 1974). Discrete, quiet, auroral arcs 
are also observed in the late evening close to the region of 
plasma flow reversal [Kamide and Rostoker, 1977]. 

Evans et al. [1977] and Maynard et al [1977] reported on 
electric fields and ionospheric currents, observed by rocket, 
near an auroral arc at the northern edge of the westward 
plasma flow. We report here on a very similar rocket flight 
over late evening auroral arcs. Arnoldy and Lewis [1977] have 
deset ibed the electron measurements on this flight and have 
shown that they can be used to predict the auroral light in- 
tensity. In this report we will describe measurements of elec- 
tric fields, ionospheric currents, and field-aligned currents in 
and near the arcs and will discuss the spatial relations between 
the arcs, the electric field, and the currents. 
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Instrumentation 

The payload included instruments to measure the electric 
field and the magnetic field (in order to detect electrical cur- 
rents) as well ns the electron flux and the auroral light. Elec- 
tric double probes were extended, on booms perpendicular to 
the spin axis, to a separation distance of 3 m. As the booms ro- 
tated in the electric field plus V X /J from rocket motion, a si- 
nusoidal signal was generated which yielded the field ampli- 
tude and direction, The electric field perpendicular to B was 
obt ained, nssuming E parallel to B to be zero, after removing 
the v X B component. The estimated accuracy of an individ- 
ual measurement of E was ±5 mV/m and ±5°, for a north- 
ward field of 50 mV/m. The sensitivity to changes in this field 
was 2 mV/m and 2°, 

Changes in magnetic field, observed as the rocket passed 
through electrical current sheets, allowed detection .and mea- 
surement of these sheets. A proton precession magnetometer, 
which provides total field magnitude measurements, wns 
flown to detect east-west horizontal ionospheric current sheets 
(electrojet currents). These sheets produce north-south mag- 
netic fields H of order 100 nT, Such fields have components H 
cos I parallel to the earth’s field (55,000 nT), where /, the field 
inclination, is 77° at Poker Flat. Changes in the parallel com- 
ponents can be observed as changes in total field magnitude. 
The proton magnetometer had an effective sensitivity, limited 
by rocket magnetic field and by signal noise, of about 10 nT. 
North-south ionospheric current sheets or field-aligned cur- 
rent sheets produce east-west megnstjo fields D, perpendicular 
to the earth’s field and therefore resulting in small magnitude 
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Fig. I. College nutgnetogram, 0600-1030 UT, Match 18, 1975, 
The period of (light 18:165 is denoted by the bar under the II trace, 
The approximate quiet time base levels of//, D, nnd Z arc shown by 
dashed lines. 


changes, ~0, 1 nT, not observable witlt the proton magnetom- 
eter. 

In order to observe more easily the relatively small changes 
in magnitude of order 10 nT, due to passage through current 
sheets, in the presence of much larger changes in magnitude 
of order 1000 nT, due to the inverse cube decrease of the 
earth’s field with altitude, it was necessary to calculate the 
earth’s field magnitude along the trajectory and to subtract 
this from the observations. Possible errors in the trajectory 
could produce a cumulative error in the difference magnitude 
MI ns large as 100 nT on the downleg. Of course, it was the 
change in M) on pnssing through the relatively thin elcetrojet 
current sheet that was essential to the current measurement 
and not the absolute accuracy of Afl. 

A sensitive flux gate magnetometer was mounted per- 
pendicular to the spin axis to detect changes in the east mag- 
netic field component due to field-aligned current sheets 
[Kin liter el ai, 1974], The times, twice per spin, when the 
transverse flux gate magnetometer (TFM) was perpendicular 
to the spin plane field component, ~25,000 nT, were mea- 
sured with an uncertainty of ±0.1 ms, Since the spin rate 
above 120 km was nearly constant, the perpendicular times 
could also be predicted. Penetration of a field-aligned current 
sheet, with resulting changes in the east magnetic field com- 
ponent, would change the direction of the spin plane magnetic 
field and thus cause a cumulative deviation of the observed 
perpendicular times from those predicted. A slow decrease in 
the spin rate during the flight, slow changes in direction of the 
earth’s field along the trajectory, and periodic variations due 
to rocket precession were modeled, and the record was cor- 
rected for these effects. After smoothing to decrease random 


fluctuations we estimated that an east component change 
greater than 200 nT, associated with a current sheet less than 
100 km thick, could be detected, Smaller or more gradual 
changes could not be clearly identified, since they might be 
due to uncertainties in modeling, The component of 8 along 
the rocket spin axis was also measured, These measurements 
agreed with the proton magnetometer results, considering the 
angle between the spin axis and the magnetic field. They were 
useful in determining the spin axis direction but did not con- 
tribute to the determination of the current sheets, The axial 
measurements will therefore not be discussed in this report. 

The electron (lux nnd auroral light measurements have 
been described by Arnold? and Lewis [1977], Electrons be- 
tween 25 eV and 15 keV were meusured with energy resolu- 
tion of 10% and flux resolution of IQ 5 (cm 3 s sr keV) 1 , The 
lower energy limit to the electron measurement is difficult to 
specify, since it depends on the unmeasured rocket potential, 
The 25-cV estimate is conservatively high, Auroral light mea- 
surements were estimated to be accurate to within 1 kR for 
the data shown. 

Launch Conditions 

Flight 18:165 was launched from the Poker Flat Rocket 
Range, Alaska, in local evening, on March 18, 1975. There 
was a positive magnetic bay (90 nT at launch) and two quiet 
auroral arcs to the north (Figures 1 and 2). The flight lasted 
over 440 s (0827-0834 UT), traveled 270 km north on mag- 
netic azimuth 4° west, and reached 224-km peak altitude, The 
spin rale after 100 s was 4.75 rps; there was a 49-s precession 
with half angle 2,4°; and the angular momentum vector was 
approximately in the magnetic meridian, 16° north of vertical 
and 29° from the magnetic field. The magnetic local time at 
launch was ~2100, 

Observations 

The auroral arcs shown in Figure 2 were clearly visible in 
the Ester Dome and Fort Yukon meridian scanning photome- 
ter records. In addition, auroral photometers mounted on the 
rocket allowed measurement of auroral light, looking directly 
down the magnetic field line below the rocket, once during 
each spin period (Figure 3). The energy flux of precipitated 
electrons was 10 ergs ern" 3 s’ 1 over the first arc nnd reached 20 



Fig. 2, The 18:165 trajectory viewed from above, Flight times are 
labeled along the trajectory. The shaded regions indicate the locations 
of auroral arcs, intense electron precipitation, and electric field ampli- 
tude decrease. The crosses show two locations of the radar drift veloc- 
ity measurements, A third location not marked is just north of Fort 
Yukon. 
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Pig. 3. Particle energy flux and auroral light, (top) The energy 
dux of precipitating electrons (a few electron volts to 15 keV) is shown 
together with the energy of a high energy peak in the electron spec- 
trum. (middle) The auroral light intensity as seen by on board pho- 
tometers looking directly down the magnetic field line is shown, (bot- 
tom) The auroral light intensity on the rocket field lines as observed 
by meridian scanning photometers at Port Yukon and lister Dome. 
This figure is a corrected version of Pigure 3 of Arnatdy and Lewis 
1 1977J, where the measured energy flux over the second arc is in error. 

over the second, while the spectral peak in the electron energy 
spectrum was less than 5 keV for the first arc and over 10 keV 
for the second. The two arcs were clearly associated with ‘in- 
verted V’ structures. The spectral peak moved beyond the en- 
ergy range of the detectors for the second arc. Calculations us- 
ing the measured spectra during this time will therefore 
represent lower limits. 

In Figure 4 the electric field is compared with the electron 
energy flux, showing that the field magnitude dropped 
abruptly as the rocket passed over the arcs where the energy 
flux was greatest. At the beginning of the flight the field mag- 
nitude was relatively high, over 60 mV/m, and pointed ap- 
proximately toward magnetic north. As the magnitude 
dropped over the first auroral arc, the direction rotated toward 
the northwest. Between the arcs the field recovered to 50 mV/ 
m. The field magnitude remained low ( — 5—15 mV/m) as the 
rocket passed over the intense arc, with direction generally in 
the southeast. After 310 s the magnitude increased slowly to 
20 mV/m, and the direction returned to northwest. Observa- 
tions of plasma flow velocity with the Chatanika radar during 
the flight allowed an independent determination of the elec- 
tric field. The radar was pointed, in sequence, toward the 
magnetic north, northwest, and northeast. Measurement inter- 
vals, along each beam direction, were chosen to provide line- 
of-sight plasma drift velocity components at 190-km altitude 
at three points along a line 170 km north of Chatanika (R. 
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pig, 4. The electric field magnitude, north component and east 
component, all in millivolts per meter, measured with the double 
probe electric field meter on 18:165. The v X B field has been sub- 
tracted. The abrupt decreases in E at 130 and 200 s occur as the rocket 
enters the regions of intense electron energy flux. Precipitated energy 
flux is shown as a histogram. 

Vondrak, private communication. 1975), The locations, 
shown in Figure 2, were just north of the bright arc, The ve- 
locity components, measured in the interval 0822 to 0833 UT, 
were used to determine the electric field, 18 mV/m north and 
4 mV/m west, The rocket measurement, north of the bright 
arc (350-390 s) was 1 5 mV/m north and 8 mV/m west. The 
agreement between the two measurements was good, consid- 
ering that the radar measurement was an average over space 
and time. 

The College magnetogram (Figure 1) indicated an eastward 
electrojet current with the current center north of College. 
The Fort Yukon magnetogram for the same period showed 
that the electrojet center was to the south of that station. Fig- 
ure 5 shows the magnetic field magnitude, measured with a 
proton magnetometer during the flight, after subtracting the 
internal reference field. The 160-nT decrease from 70 to 110 s 
was caused by penetration of the eastward electrojet current 
layer, while the relatively slow return of the magnitude to aero 
from 300 to 400 s suggested that the electrojet was not directly 
penetrated on the downward portion of the flight. The magni- 
tude change, 160 nT, corresponds to a horizontal component 
H change of 700 nT, 



O IOO 200 300 3 00 


FU0HT TIME (sue) 

Fig. 5. Magnetic signature of the eastward electrojet. The magni- 
tude of the earth’s internal field has been subtracted from the proton 
precess'on magnetometer record to emphasize the field changes due to 
ionospheric currents. The abrupt change between 70 and 1 10 s in- 
dicates passage upward through an eastward directed electrical cur- 
rcntshecL The absence of a similar abrupt transition in the downward 
flight record suggests that the rocket came down to the north of the 
electrojct northern boundary. 
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In addition to the eastward electrojet current, field-aligned 
currents carried by electrons (a few electron volts to 15 keV) 
were measured, as shown in Figure 6. The net field-aligned 
current was calculated from the measured upward and down- 
ward electron fiuxes. The upward (conventional) current den- 
sity over the weaker arc, 4 X 10 6 A/in 3 , was about twice that 
over the more intense arc, 2 X liT 6 A/m 3 , The average cur- 
rent over the 100-km-thick current sheet was 2.5 X IQ" 6 A/m 3 , 
This should produce a total transverse (east to west) field 
change of 320 nT, The direct particle measurements of field- 
aligned currents were supplemented by sensitive measurement 
of the direction of the magnetic field component per- 
pendicular to the rocket spin axis by the transverse flux gute 
mugnclometcr. The cumulative change, between 150 and 300 
s flight time, in the time, each spin, Of observing; the transverse 
magnetic field was 2,6 X 10” 3 s and was equivalent to a 
change (east to west) in magnetic field of over 1000 nT, as 
shown in Figure 6, Assuming this field change resulted from 
penetration of an upward flowing field-aligned current sheet 
(extending east-west) 100 km thick, the (uniform) current den- 
sity was ~8 X 10 ‘ A/m 3 : 

Jw “ ■— A/m 3 
floA y 

where Afl c „„ is the change in the east component shown in 
Figure 6 and Aj' is the sheet thickness, This measurement of 
field-aligned current gives a result about 3 times greater than 
the particle measurements. The disagreement between the two 
techniques for measuring field-aligned current could be due to 
errors in the modeling or to field-aligned current carried by 
ions or thermal electrons not measured by the particle detec- 
tors. 

C \CULAT10NS 

Height-integrated Hal: ,nd Pedersen conductivities arc 
shown in Figure 7, These were calculated using the incoming 
measured electron flux, the neutral model atmosphere of Rees 
el ah [1977], the temperature profile of Jones and Rees [1973], 
and the ion composition of an auroral arc given by Rees and 



Fig, 6, Field-aligned current from electron flux and magnetic 
measurements. The difference between upward and downward mov- 
ing electrons is a net downward electron flux. This constitutes the up- 
ward electrical current density shown in the bottom panel, The flight 
photometer record (5577 A) is also plotted to show the relation of 
field-aligned current to the auroral arcs. In the top panel the change 
in the eastward magnetic field component, as deduced from the trans- 
verse flux gate measurements, is shown. 



Fig. 7. Height-integrated Hall and Pedersen conductivities calcu- 
lated from the precipitating electron flux, 


Walker [1968], First the electron production rate r/(Z) result- 
ing from the electron flux impinging On the utmosphere at 
each altitude level was calculated according to the work of 
Rees [1963], Equation 1 in Rees’s paper was modified to apply 
to a measured directional intensity of electrons J(E), as fol- 
lows: 

q(Z) ~2tt J J J(E) sin 0 cast) ME 

where the ratio qz/F is equation 1 of Rees [ 1963], using the 
isotropic, normalized energy dissipation distribution function 
of that paper, The electron density was obtained from the pro- 
duction rate under equilibrium conditions by using the effec- 
tive recombination coefficients tabulated by Evans et ah 
[1977], Finally, the collision frequencies and the Hall and 
Pedersen conductivities were obtained as outlined by Evans et 
ah [1977], Note in Figure 7 that the Pedersen conductivity £ ; , 
exceeds the Hall conductivity 2 W for the first arc. The reverse 
is true for the second and brighter arc, which is produced by 
more energetic electrons that deposit their energy deeper in 
the atmosphere, 

One may calculate the Joule power dissipated by the iono- 
spheric current component parallel to E from the calculated 
Pedersen conductivity and the measured electric field as 2,, E 2 
[Evans el ah, 1977]. This Joule power as a function of flight 
time is given in Figure 8 along with the power carried into the 
atmosphere by the precipitating electrons, The Joule power 
decreased abruptly as the rocket moved over the aurora due to 
the decrease in the electric field, while the precipitated power 
peaked over the arcs, The sum of the two was about the same 
over the arcs as between or south of them. Note that very little 
power was present north of the auroral arcs, suggesting that a 
natural boundary was crossed. 

The horizontal electrical current can also be obtained by 
using the calculated conductivities and the measured electric 
field as of nonh ~ — 2//E eus , and */„„„ , — 4* 

2 // £'„ 0 „ h . The calculated current components are shown in 
Figure 9. In the lower section the eastward current density is 
substantial, ~0,6 A/m from 90 to 110 s. It then decreases 
steadily to ~0.2 A/m at 140 s. This calculated component is 
weak and erratic in the first auroral arc, recovers to 0.2 A/m 
between arcs, and is weak and occasionally westward in the 
second arc. The northward component in the upper section 
continues at —0,4 A/m from 90 to 200 s through the first arc, 
then decreases to —0.4 A/m at 280 s. It increases rapidly to 0.2 
A/m at 310 s. The uncertainty in these calculated current den- 
sities is principally due to the uncertainties in the ionosphere 
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Fig, 8, Powet dissipation in the ionosphere. The energy per sec- 
ond deposited by the precipitating electrons is highest in the auroral 
arcs. The Joule power is strongest outside the nrcs where the electric 
field is higher. The sum of these two decreases abruptly north of the 
bright arc (300 s). 


modeling. We estimate that the current densities may be in er- 
ror by as much ns ±50%, 

Discussion 

These rocket measurements provide an accurate com- 
parison between the locations of several auroral zone features: 
the northern edge of the westward plasma (low, bands of in- 
tense electron precipitation and the associated auroral arcs, 
horizontal ionospheric electric currents, and field aligned cur- 
rents. We emphasize that the auroral arcs and magnetic field 
were very stable during the flight, so we interpret observed 
changes as due to rocket passage through and over stationary 
structures. 

The two tires were produced by observed fluxes of energetic 
electrons. A moldy and Lewis [19 7 7] have shown that the au- 
roral light can be predicted from the observed electron fluxes 
and vice versa. Two field-aligned sheets of precipitating elec- 
tron flux stretched from the western to the eastern horizon 
with bands of auroral light at the lower edges, 

The electric field observations show that the auroral arcs 
were imbedded in the northern edge of a strong westward 
plasma flow. The decreases in the electric field over the au- 
roral arcs imply polarization charge layers at the edges of the 
arcs. Similar decreases in electric field over auroral arcs were 
seen in earlier flights [Potter, 1970; Maynard et al, 1973). 
Sometimes the electric field docs not decrease over arcs, as, 
for example, during a flight reported by Kintner et al, [1974] 
over a rapidly moving arc in auroral breakup. It appears that 
only in stable, quiet arcs, like the evening arcs discussed here 
or the similar arc reported by Maynard et al, [1977], are the 
polarization charge layers established to reduce the electric 
field. 

We must rely on ground magnetic observations and the 
electric field/conductivity calculations as well as rocket obser- 
vations to deduce the distribution of the horizontal iono- 
spheric currents. A strong eastward current layer is detected 
by the proton magnetometer in the upward leg of the flight 
(40-50 km downrange), while the current is weak or absent on 
the downward leg (280-290-km range). This establishes that 
the eastward electrojet terminated between 50 and 280 km 
north of Poker Flat. The height-integrated current density for 


the eastward clectrojet derived from the change in field mag- 
nitude is All/fto cos / -> 0.57 A/(m in the N-S direction), in 
good agreement with the calculated value (from measured E 
and calculated conductivity) of 0.6 A/m at about the same lo- 
cation (Figure 9). Encouraged by this agreement, we note that 
the calculated eastward current density decreases to zero in 
the first arc, recovers to 0,2 A/m between arcs, and stays near 
zero thereafter. These calculated results indicated that the 
eastward clectrojet is greatly reduced in the first arc, 80 km 
north of Poker Flat, and terminated in the second arc, 175 km 
north. A third indication of the electrojet extent is obtained 
from ground magnetometers. Figure 10 shows A// and A Z for 
several observatories in the Alaskan sector; A H indicates an 
eastward clectrojet with maximum current near Poker Flat 
and College, and AZ crosses zero just north of Poker Flat, in- 
dicating the center of the electrojet current lies there. The 
shape of A// and AZ curves shows that the current is confined 
to the region near Poker Flat with the northern edge near Fort 
Yukon. The southern edge is less certain because of lack of 
observations between College and Sitka, These ground mea- 
surements don’t allow a precise determination of the northern 
edge of the eastward electrojet, but they do support the rocket 
magnetometer and the calculated results. 

Our information about northward current comes solely 
from the calculations. We had hoped the TFM measurements 
could be used for these currents, but the rocket attitude 
changes too rapidly below 120 km to allow reliable modeling. 
The D component ground magnetometer results in Figure 1 
(<25 nT at launch) indicate small and decreasing northward 
current. The calculated value stuys above 0.4 A/m, however, 
until entry into the second arc 120 km downrange. For a wide 
northward horizontal current sheet the D component below 
the current sheet is D = lkJn/1 = 250 nT, 10 times the ob- 
served D component on the ground, The discrepancy is prob- 
ably due to the solenoidal nature of the meridional currents, 
downward field-aligned currents south of Poker Flat, north- 
ward current south of the arcs, and upward field-aligned cur- 
rents over the arcs, with the magnetic field confined within the 
solenoid. 

There have been many observations of field-aligned cur- 
rents with polar-orbiting satellites and by rocket [Armstrong 
and Zmuda, 1973; Potemra, 1979; Cloutier and Anderson, 
1975]. Current densities from 1 to 10 jttA/m 2 have been re- 



Fig, 9. Horizontal current calculated from electric field and con- 
ductivity, The eastward current in amps per meter; the integrated 
current from top to bottom of the clectrojet in a t-m-thick sheet, 
aligned cast-west) is the component that can be detected in the proton 
precession magnetometer record. 
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Fig* 10. Profile of the eastward electrojet as seen in the per- 
turbations in the northward and vertical magnetic field components at 
several Alaskan observatories. Uncertainties in determining these per- 
turbations from ground magnetograms arc about ±20y. 


ported. The rocket flight discussed here provided three meth- 
ods of determining the field-aligned current. The current re- 
sulting from the net flux of 25 Ev-15 keV electrons over the 
auroral arcs (2.5 /iA/m 5 average) is shown in Figure 6 together 
with the change in cast-west magnetic field. The field change 
implies an average current of 8 pA/tn 2 , Over the ~I00-km- 
thick sheet of upward current the integrated current density is 
-0.8 A/m. This agrees with the decrease, from +0.4 to -0.4 
A/m, in the calculated northward current (Figure 9). Down- 
ward field-aligned current is also implied by several increases 
in northward current density, particularly from 280 to 310 s in 
Figure 9. Since both the magnetic measurements and the di- 
vergence of the northward ionospheric current give an Up- 
ward field-aligned current larger than the measured electron 
current (25 eV-*15 keV), we attempted to extrapolate the elec- 
tron measurements to tower and higher energies than those 
measured, Starting with the measured electron flux at 25 eV, 
the flux was extrapolated down to 0. 1 eV, using a power law 
decrease of flux with energy E~'. The flux of these low-energy 
electrons is assumed to be isotropic, equal intensities upward 
and downward. However, in order to find the maximum pos- 
sible current due to these unmeasured low energy electrons we 
neglected the upward streaming electrons, This calculation 
gave an additional upward current density of 2.5 pA/m 2 . Ex- 
trapolation upward in energy from 15 keV to 50 keV in a sim- 
ilar fashion gave less than I pA/m 2 , We estimate that the av- 
erage field-aligned current carried by precipitating electrons, 
including these extrapolated values, was less than 6 pA/m 2 , 
The remaining discrepancy is perhaps due to unmeasured up- 
streaming thermal positive ions, to uncertainties in the mag- 
netic measurements, or to uncertainties in the calculated 
northward current. 

One further source of error may be neutral winds, particu- 
larly in the interval 210-280 s, when the electric field was less 
than 10 mV/m. Neutral wind studies by Brekke el al. [1973] 
with the Chatanika radar indicate that neutral winds near 
2100 local time are southeast and weak when strong con- 
vective ion flow (strong electric field) is absent. In the present 
case there is, however, strong westward ion flow, thus west- 
ward neutral wind due to ion drag, just south of the region 
of weak electric field and weak ion flow. One might expect, 
from continuity, that the neutral winds in the weak field re- 
gion would be south, A weak southward neutral wind should 
produce a weak southward ion flow by collisions and thus a 
southward current. This additional current from the neutral 
wind would add to the calculated southward current in Figure 


9 (210-280 s) and would thus increase the upward field- 
aligned current in this interval. 

A model current system is shown in Figure 1 1 with the east- 
ward, northward, and upward field-aligned components dis- 
cussed above. Diffuse downward field-aligned currents pre- 
sumably exist south of Poker Flat, but these have not been 
included in the model. Field magnitude A# and components 
A//, A D, and AZ produced along the rocket trajectory by this 
model current arc also shown. The A B and AD agree approxi- 
mately with the proton magnetometer and transverse flux gate 
magnetometer results, The observed A B of Figure 5 docs not 
show the rapid recovery of the model from 100 to 200 s. 
Bringing the termination of the model eastward current far- 
ther north would improve the agreement; AH, AD, and AZ 
agree reasonably well at zero seconds with the Poker Flat 
ground observations. Note that AD is ~0 except above 100 km 
in the solenoid, Although the model AH at zero seconds, 100 
n'f, agrees approximately with the ground magnetometer, 90 
nT, both arc much less than the field, 700 nT/2 - 350 nT, di- 
rectly beneath the electrojet, inferred from the proton magne- 
tometer assuming an infinite current sheet. In this case, since 
the current sheet is relatively narrow, AH is only 200 nT be- 
neath the electrojet, as seen in Figure 1 1, and a change in Z as 
well as in H contributes to the observed change in B as the 
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Fig. II. Magnetic field components from model current calcu- 
lations. The model is shown in the upper left panel. The AH is the 
magnetic north component principally due to the eastward electrojet 
current, AD the magnetic east component due to northward and field- 
aligned currents, AZ the vertical component due also to the eastward 
current, and A B the change in total field due principally to the east- 
ward current, The AA, the component along the rocket spin axis, was 
used to compare with the axial fiux gate magnetometer records not 
discussed in this report. The model consists of an eastward current, 12 
pA/nt 2 , extending front 30 km south of the launch site to 70 km north 
and from 100 to 150 km in altitude; a northward current, 8 pA/nt 2 , ex- 
tending from 300 km south of launch to 100 km north and 100 to 150 
km in altitude; and a field-aligned current, 10 p A/m 2 , extending from 
75 to 150 km north of launch and 120 to 700 km altitude. All currents 
were essentially infinite east-west. 
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rocket passes through the current sheet. At 440 s, A// and AZ 
are considerably less than the Fort Yukon values in Figure 10. 
Rocket impact was about 70 km north of Fort Yukon, so they 
should be somewhat less. Again, extending the model east- 
ward current farther north would increase A// and iiZ at 440 
s. Although the model could be improved, it is adequate to 
show the essential agreement between model, rocket, and 
ground observations. For example, addition of an induced 
subsurface current parallel to the clcctrojet would increase A If 
and decrease AZ of the model at ground level by as much as 
50CF but would have less effect on the fields along the upper 
rocket trajectory. 

There is a substantial energy deposit in the auroral iono- 
sphere over the segment traversed by the rocket. The power 
level by precipitating electrons over the arcs and by Joule 
heating of the northward current south of the arcs remains 
near 15 ergs/cnr s until the northern edge of the second au- 
rora. This energy flux (15 ergs/enF s ■= 0.015 W/nr) amounts 
to 1.5 MW for the segment, 1 km east-west and 100 km north- 
south, that the rocket (lew over 

Figure 12 summarizes the findings of this flight. The rocket 
entered the ionosphere at 60 s in a strong westward plasma 
(low and eastward electrojet (Hall current). The northward 
electric field, westward plasma flow, and eastward Hall cur- 
rent weakened over one auroral arc and terminated over a 
second, more intense auroral arc. Field-aligned currents were 
observed streaming out of the auroral arcs, providing continu- 
ity with northward ionospheric currents. The second arc ap- 
peared to be a very definite boundary to the plasma llow/cur- 
rent region in that all observed parameters were weak and 
erratic north of the arc. We are tempted to identify this region 
as the Hnrang discontinuity [Maynard, 1974; Maynard el al., 
1977], There is no evidence for a westward electrojet, how- 
ever, and the observations are, at 2100 MLT, a little far west 
for the Hnrang pattern. There is no clear evidence of a large- 
scale eastward plasma flow (southward electric field) in the 
rocket observations so it isn’t clear the rocket entered the po- 
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Fig. 12. Schematic view of the electric field vectors al 20-s inter- 
vals along the trajectory. Ex U drift velocity (obtained by multi- 
plying E by 20 and rotating counterclockwise 90°) and calculated hor- 
izontal current are also shown. The scale at right shows the 
approximate horizontal range of the rocket. 


Inr cap. The region north of the arc may be a local eddy in the 
plasma flow. 

The observations provide a strong confirmation of the ear- 
lier similar observations of Maynard et al. [1977] and Evans ct 
al. [1977], These two sets of observations provide a detailed 
description of the northern edge of the westward plasma (low/ 
eastward electrojet/quiet arc region in the evening hours, 

Acknowledgments. We acknowledge the assistance of Steve Wer- 
kema, John Lebens, and Wai Kwok Kwong with data processing and 
of John Humenansky with payload design and construction at the 
University of Minnesota. At the University of New Hampshire the ef- 
forts of Mark Widholm. Peter Lewis, and Ralph Varney are appreci- 
ated. The University of Minnesota was supported by NASA grants 
NGR24-005-253 and NSG-6G05, and the University of New Hamp- 
shire by grants NGR30-002-054 and NSG-6013. A portion of this 
work at TRW was supported by NASA contract NASW-3087. We 
thank the College Observatory and World Data Center A for ground 
magnetograms. Laurence Cahill is grateful for the support of Gerhard 
Haerendcl and the pec le at the Max Planck Institute for Extra- 
terrestrial Physics where this report was completed. Finally, we wish 
to acknowledge the support given by Richard Vondrak and the Chat- 
anika radar stall’ for this rocket flight. The Chatanika Incoherent Ra- 
dar Facility is operated by the Stunford Research Institute and sup- 
ported by the Defense Nuclear Agency (control number DNA001- 
74C-QI67) and the National Science Foundation (grant DES72- 
01644-A02). 

The Editor thanks II. Anderson and M. H. Rees for their assistance 
in evaluating this paper. 


References 

Armstrong, J. C., and A. J. Zmtida, Triaxial magnetic measurements 
of field-aligned currents at BOO km in the auroral region: Initial re- 
sults, J. Geophys. Res., 78, 6802, 1973, 

Arnoldy, R. L.. and l 1 . B. Lewis, Jr„ Correlation of ground-based and 
topside photometric observations with auroral electron spectra 
measurements at rocket altitudes, J. Geophys. Res., 82, 5563, 1977. 

Brekke, A., J. R. Doupnik, and P. M. Banks, A preliminary study of 
the neutral wind in the auroral E region, J. Geophys. Res., 78, 8235, 
1973. 

Cloutier, P. A., and II. R. Anderson, Observations of Birkeland cur- 
rents, Space Sci. Rev., 17, 503, 1975. 

Evans, D. S., N, C. Maynard, J. Trpitn, T. Jacobsen, and A. Egeland, 
Auroral vector electric field and particle comparisons, 2, Elec- 
trodynamics of an arc, J. Geophys, Res., 82, 2235, 1977. 

Jones, R. A., and M. II. Rees, Time-dependent studies of the aurora, 
1, Ion density and composition, Planet, Space Scl„ 21, 537, 1973. 

Knmide, Y, F„ and G. Rostoker, The spatial relationship of field- 
aligned currents and auroral electrojets to the distribution of night- 
side auroras, J. Geophys. Res., 82, 5589, 1977. 

Kintner, P. M., L. J. Cahill, Jr., and R. L. Arnoldy, Current system in 
an auroral substorm, J. Geophys. Res., 79, 4326, 1974. 

Maynard, N. C., Electric field measurements across the Hnrang dis- 
continuity, J. Geophys. Res., 79, 4620, 1974. 

Maynard, N. C., A. B-hnsen, P. Christophcrson, A. Egeland, and R. 
Lundin, An example of anticorrelation of auroral particles and 
electric fields, J, Geophys, Res., 78, 3976, 1973. 

Maynard, N. C.. D. S. Evans, B. Machlum, and A. Egeland, Auroral 
vector electric field and particle comparisons, 1, Premidnight con- 
vection topology, J. Geophys. Res., 82, 2227, 1977. 

Potemru, T. A., Current systems in the earth’s magnetosphere, Rev. 
Geophys. Space Phys., 17, 640, 1979. 

Potter, W. E., Rocket measurements of auroral electric and magnetic 
fields, J. Geophys. Res., 7S, 5415, 1970. 

Rees, M. H., Auroral ionization and excitation by incident energetic 
electrons, Planet. Space Sci., 11, 1209, 1963. 

Rees, M. II., and J. C. G. Walker, Ion and electron heating by auroral 
electric fields, /1/m. Geophys,, 2d, 193, 1968. 

Rees, M. H., A. I. Stewart, W. E. Sharp, P. B. Hays, R. A. Hoffman, 
L. H. Brace, J. P. Doering, and W. K. Peterson, Coordinated rocket 
and satellite measurements of an auroral event, 1, The satellite ob- 
servations and analysis, J. Geophys. Res., 82, 2250, 197T 

(Received June 1, 1978; 
revised March 17, 1980; 
accepted March 18, 1980.) 





»'0L. 84, NO. A9 


JOURNAL OF GEOPHYSICAL RESEARCH 


SEPTEMBER 1, 1979 


IONOSPHERIC ELECTRICAL CURRENTS IN THE LATE EVENING PLASMA FLOW REVERSAL 

D. A. Behm, 1 F. Frimdahl, 2 L. J. Zanetti, Jr., 3 ’ 1 ' 

R. L. Arnoldy, 3 and L. J. Cahill, Jr. 1 


Abstract . An instrumented sounding rocket was 
launched from Anddya, Norway, in January 1977, in 
the late evening auroral oval. Its trajectory 
took it northward over a quiet auroral arc and 
over a region of plasma flow reversal. Horizon- 
tal currents were inferred from magnetic field 
measurements; they were also calculated from the 
measured electric field and conductivities deriv- 
ed from energetic particle influxes. South of 
the reversal region the sheet current density of 
the eastward electrojet was found to be '\-0.1 A/m, 
while north of the reversal the current density 
was less than 0.05 A/m. The currents were enhanc- 
ed in the arc. Current continuity considerations 
in the meridian plane implied downward field- 
aligned currents at the southern edge of the arc 
and upward field-aligned currents at the northern 
edge. 

Introduction 

In January 1977 a Bounding rocket, launched 
from Anddya, Norway, (69.3°N, 16.0°E geographic; 
66.5°N invariant magnetic latitude) traversed a 
region, in the late evening auroral oval, where 
the plasma flow reversed from westward to east- 
ward. Such plasma flow reversal, often called 
the Harang discontinuity, has been observed pre- 
viously through barium ion drift, electric field 
observations, and radar Doppler measurements 
[Wescott et al., 1969; Maynard, 1974; Wedde et 
al., 1977]. In earlier reports on this rocket 
flight we described the precipitating electron 
flux and its relation to an auroral arc in the 
southern part of the flow reversal [Arnoldy, 

1977] and the flow reversal itself as simulta- 
neously observed by electric field double probes 
on the rocket and by the Stare (Scandinavian Twin 
Auroral Radar Experiment) radar [Cahill et al . , 
1978b] . In this report we describe horizontal 
ionospheric currents, flowing near the flow re- 
versal, as inferred from ground and rocket mag- 
netic observations and as calculated from derived 
ionospheric conductivities and the observed elec- 
tric field. 

Instruments and Trajectory 

The instrumentation, which Included several 
magnetometers, electric field double probes, and 
curved plate electrostatic analyzers, has been 
described in earlier reports [Cahill et al., 
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1978b; Arnoldy and Lewis, 1977]. Two magnetom- 
eters were flown to detect the east-west compo- 
nent of the horizontal ionospheric currents. A 
proton precession magnetometer (PPM) provided 
total field magnitude measurements once each 
second with an accuracy, limited by vehicle 
noise, of 'VlO nanotesla (nT) . A three-component 
flux gate magnetometer was also flown to supple- 
ment the PPM measurements and to evaluate the 
flux gate as a possible replacement for the PPM 
on future flights. Three perpendicular component 
measurements were sampled about 360 times each 
second with a resolution of 30 nT. Only the 
field magnitude, obtained by combining the three 
components, was used in this report. Both in- 
struments were sensitive, as described below, to 
changes in the north-south and vertical field 
components observed as the rocket penetrated 
east-west current layers. 

The magnetic fields produced by the horizontal 
current layers were of order 100 nT and horizon- 
tal (except near the layer edges) , while the 
earth's field at Anddya was about 55,000 nT and 
only 12° from vertical. Magnetic north was 10° 
west of geographic north. An eastward (or west- 
ward) flowing horizontal current layer produced a 
northward (or southward) magnetic field (X) di- 
rectly below the layer. Such a field made a 
contribution to the total field magnitude of 
X, sin 12° (20.8 nT when X was 100 nT) . This con- 
tribution changed sign as the rocket climbed to 
the topside of the current layer. A northward or 
southward flowing current layer produced an east- 
ward or westward magnetic field also of order 100 
nT but perpendicular to the earth's field. The 
total field magnitude increased in this case by 
only about 0.1 nT, an undetectable change. 

The electric field double probes provided 
measurements of the north and east electric field 
components with an accuracy of +5 mV/m [Cahill et 
al., 1978a]. Electron flux measurements, between 
a few electron volts and 22 keV, were available 
with energy resolution of 10% and a flux resolu- 
tion of 10“ (cm 2 s sr keV) -1 . With the use of 
multiple detectors, distribution functions were 
measured every 5 s of flight time [Arnoldy and 
Lewis, 1977]. 

The rocket was launched at 1918:06 UT ('«2200 
MLT) on January 23, 1977. The flight lasted 470 
s, reached 221-km peak altitude, and travelled 
213 km downrange on geographic azimuth 345°. The 
spin rate was 3,9 Hz, and there was a precession 
of 3.24° half angle with period 58.55 s. 

Observations 

A positive magnetic bay (40 nT) in the north- 
ward geographic component, an upward vertical 
component (20 nT) , and almost zero disturbance in 
the eastward geographic component existed at 
launch, as shown in Figure 1, while a quiet auro- 
ral arc was visible in the north. Observers at 
the Stare radar reported a westward plasma flow 
over Anddya and an eastward flow about 100 km 
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Fig. 1, Anddya magnetogram, X is the north, 

Y the east, and Z the vertical (down) compo- 
nent, The vertical lines indicate launch and 
reentry times. The dashed portion of each 
trace represents a time period when the trace 
was lost owing to interference. 

further north; the conditions were right for 
flight over the Harang discontinuity [Cahill et 
al., 1973b] . 

j The magnetic signature of penetration through 
an eastward electrojet is shown in Figure 2. The 
proton precession magnetometer record is shown on 
top and the flux gate record at the bottom. Sim- 
ilar reference fields have been subtracted in 
each case to emphasize the changes due to pene- 
trating the currant layers. Between 90 and 120 
km on the upleg a sudden 50-nT decrease in AB was 
apparent in both records. There were several 
other smaller features (a.lO nT) in each record, 
presumably not due to ambient magnetic field 
changes, since there were no correlations between 
the two records. Between 120 and 90 km on the 
downleg there was no field decrease that could be 
attributed to a possible westward electrojet cur- 
rent layer. The 50-nT change in AB corresponded 
to a 50/ sin 12° » 240 nT change in the X compo- 
nent, For an infinite current sheet, one half of 
this, 120 nT, would be observed on the ground 
below. The AX » 240 nT change in the horizontal 
component through the current layer corresponded 
to a sheet current strength of AX/po , v 0.19 A/m, 
The electric field measurements, shown as the 
top panel of Figure 3, have been reported by 
Cahill et al. [1978b]; a brief review follows. 
South of the auroral arc a strong northward elec- 
tric field produced a westward plasma flow and an 
eastward (Hall current) electrojet. The field 
decreased over the arc, decreased again north of 
the arc, and rotated through an eastward direc- 
tion (northward plasma flow) . Recovering north 
of this, the field rotated to the south (eastward 
plasma flow) . A sudden increase to over 50 mV/m 
at 300 s was followed by a direction shift to 
southeast, an oscillation with 6-s period, and a 
gradual decrease in magnitude. An eastward elec- 
tric field (northward plasma flow) is not nor- 
mally associated with plasma flow reversal in the 
Harang discontinuity, It is possible that the 
weak eastward field is a transient feature asso- 
ciated with counterclockwise flow around an auro- 


ral loop observed east of the trajectory [Cahill 
et al., 1978b]. 

Electron measurements on this flight have also 
been reported earlier [Arnoldy, 1977]. Briefly, 
as is shown in^Figure 4, there were intense, i0* 
(cm 2 s sr keV) 1 " 1 , electron fluxes with energies 
of a few hundred electron volts over and near the 
arc (130-160 s) , in the plasma flow reversal re- 
gion (200-210 s) and in the eastward plasma flow 
region (250-280 s). Electrons of keV energies 
were observed only directly over the arc. The 
electric field magnitude, inverted, is also 
shown. Note that the field ip low when the 
0.3-keV flux is high. 

Calculations 

The east Jg and north components of the 
horizontal vertically integrated current density 
were calculated [Anderson and Vondrak, 1975; 

Evans et al., 1977; Cahill et al., 1978a] . The 
electron flux and energy spectrum were used, with 
models of the ionosphere, to calculate electron 
densities [Rees, 1963]. Height-integrated Hall 
and Pedersen conductivities, shown in Figure 3, 
were obtained from the electron densities and 
model collision frequencies. In turn, the con- 
ductivities and measured north and east electric 
field components Ejj and Ee yielded north and east 
ionospheric current components Jn and Je* 

The height-integrated conductivities were 
about 2 siemens (S) (1 S 5 1 mho) ir. the westward 
plasma flow region, south of the auroral arc. The 
Hall conductivity increased to 20 S over the arc, 
while the Pedersen conductivity went to 32 S. 
North of the arc the calculated conductivities 
remained below 2 S. Since the northward electric 
field was high in the strong westward plasma flow 
south of the arc, there was a moderate calculated 
current density there, , v0.1 A/m for the eastward 
component and ^0.05 A/m for the northward compo- 
nent. The current densities increased dramati- 
cally in the auroral arc, to a.o.S A/m east and 
, v0.4 A/m north. North of the arc, including the 
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Fig. 2. Magnetic field magnitude measure- 
ments . Reference fields have been subtracted 
from the measured fields. The zero levels are 
arbitary. 
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strong eastward plasma flow from 240 to 300 s, 
the calculated current density components were 
less than 0.05 A/m. 

The Joule power dissipation by horizontal ion- 
ospheric current flowing in the direction of E 
was also calculated [Evans et al., 1977]. The 
Joule power density is compared, in Figure 3, 
with the power deposited in the same column of 
the ionosphere by precipitating electrons. The 
Joule power was a> 3 mW/m 2 column (1 mW/m 2 « 1 
erg/cm 2 s) in the eastward qlectrojet south of 
the arc, principally due to the northward current 
component. Over the arc the Joule dissipation 
increased to <v8 mW/m 2 column, while north of the 
arc the power was less than 1 mW/m 2 column. The 
precipitating electron power was less than 1 
mW/m 2 , except over the arc, where it rose to 7 
mN/m 2 . 

Discussion 

We have two measurements of the magnetic field 
due to the eastward electrojet current. The 
northward component at Anddya is AO nT and the 
northward component just below the current layer 
is 120 nT. For a very broad current sheet Che 
full 120 nT should be seen at ground level. The 
section of the eastward electrojet penetrated by 
the rocket must therefore be of limited width 
from north to south. At Anddya the vertical com- 
ponent is -20 nT, indicating that the center of 
the eastward electrojet is south of Anddya. 

An east-west band current of 0,19-A/ra strength 
(Inferred from the rocket measurements) approxi- 
mately 150 km south of the launch site would ex- 
plain the 40-nT horizontal disturbance and the 
-20 nT vertical disturbance observed at the 
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Fig. 3. Measured electric field, calculated 
conductivities, calculated horizontal currents, 
and power dissipated by electrons and by Joule 
heating. The electric field and horizontal 
current components are given with reference to 
magnetic north and east. The arc is located 
between 137 and 152 s. 



Fig A. Electron directional intensities at 
0.3 and 3.0 keV. The higher energy electrons 
are, seen only over the arc (137 to 152 s) . 

The electric field magnitude is shown with 
Inverted scale at right. 

ground. The height-integrated current density 
derived from the rocket magnetometer measurement 
in 0.19 A/m, while the eastward current density 
calculated in the eastward electrojet, 20 km 
south of the arc but 25 km north of the rocket 
electrojet penetration, is *'■0.1 A/m. The calcu- 
lated current is about half that derived from the 
magnetometer measurements. This may be due to 
variations in current density at the two differ- 
ent locations but may also be caused by uncer- 
tainties in the modeling procedures used to cal- 
culate the current. Considering the uncertainties 
in various models used, we estimate the possible 
error in the conductivities to be less than +50%. 
Since the errors in the electric field are 
smaller, the calculated currents may also be in 
error by as much as +50%. 

The currents were"~calculated on the assumption 
that the neutral winds were negligible. No neu- 
tral wind measurements were available during the 
flight; however, observed neutral winds in the 
auroral zone during geomagnetically quiet or mod- 
erately disturb ad times (Kp <_ 3) have been found 
by Brekke et al. [1973] to be generally less than 
or of the order of 100 m/s. Neutral winds of 
this magnitude, if they were present, would add 
an effective electric field vector to the observ- 
ed electric field vector of no more than 5 mV/m. 
The effects of neutral winds on the calculated 
currents would be small except in the reversal 
region, where the measured electric field is 
small. 

The lack of a westward electrojet as observed 
by the magnetometers on the downward leg agrees 
with the calculations. A very weak, 0.02-0.0A 
A/m, calculated westward component is present 
below the rocket from 250 to 300 s. With the 
electric field direction change to southeast the 
calculated currents become southward, 0.02 to 
0.0A A/tn, 

The electric fields, corresponding plasma flow 
velocities, and calculated horizontal current 
vectors are shown in a summary view in Figure 5. 
The weakening of the electric field and of the 
plasma velocity over the arc is apparent as well 
as the reversal of the flow velocity north of the 
arc. The largest current densities, over 0.5 A/m 
to the northeast, are associated with the arc. 
Unlike earlier flights over quiet arcs, in this 
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Fig. 5. Summary plot of electric field, 
plasma convection velocity, and calculated 
horizontal currents. The shaded region de- 
notes the location of the arc. The trajectory 
has been projected to 110 km; therefore the 
time scale on the left is not linear. Scales 
for the three quantities are given at the 
bottom. 


case the calculated eastward electrojet current 
is intensified in the arc [Evans et al., 1977; 
Cahill et al., 1978a]. This is principally be- 
cause the conductivities increase by & factor of 
10 over this arc but only by a factor of 3 or 4 
over the earlier arcs. 

The total current in the arc is about 4 kA, 
and the payload flies by it at a distance of 
about 25 km. The peak magnetic disturbance in 
the total field at the rocket is then expected to 
be of the order of 30 nT, which should be detect- 
able by the on-boatd magnetometers. However, the 
disturbance increases and decreases gradually 
over a substantial part of the trajectory, so 
with the modeling technique used in reducing the 
magnetometer data the disturbance may not be dis- 
cernable. The lack of evidence of the line cur- 
rent to the magnetometer measurements could, of 
course, also be due to errors in the ionospheric 
modeling procedure lending to a calculated arc 
current that is too large. As was estimated ear- 
lier, the errors may be as large as +50%. 

The electric field was almost constant during 
the first part of the flight from 83 to 116 s, 
directed ''>6° west of geographic north. During 
this time interval the And(iya horizontal magnetic 
perturbation was also constant and directed to- 
ward M* west of geographic north or, within the 
measurement uncertainties, parallel to the dc 
electric field. With calculated height-integra- 


ted Hall and Pedersen conductivities of S in 
this region the northward Pedersen current should 
have rotated the horizontal magnetic disturbance 
45* to the west. The missing Pedersen current 
magnetic field component at ground level confirms 
Vasyllunns' [1970] suggestion that the magnetic 
fields from the field-aligned currents and the 
Pedersen currents are confined in regions above 
the ionosphere. The northward current must be 
part of a solenoids! current system where current 
flows down in a field-aligned east-west sheet, 
then north horizontally, then up in a field- 
aligned sheet, 

Homogeneity along the auroral arc (east-west) 
being assumed, changes in the northward current 
density imply field-aligned currents in order to 
satisfy current continuity [Evans et al,, 1977; 
Cahill et al, , 1978a]. The Increase from 0.03 
A/m at 90 s to 1.0 A/m at 120 s Indicates a down- 
ward field-aligned current sheet 20-km thick, 
with an integrated current density of 0.07 A/m, 

In similar fashion a 0,2-A/m downward sheet is 
implied in the southern part of the arc, and a 
0.3-A/m upward Bhect in the northern part. 

We are completing a study of field-aligned 
currents measured by several techniques during 
this flight: net electron fluxes (from a few cV 

to 22 keV) ; ion drift measurement and measure- 
ments of changes in the east-west magnetic compo- 
nent when passing through possible field-aligned 
current sheets, There appear to be some discrep- 
ancies between the field-aligned currents meas- 
ured by the several techniques and between the 
measured currents and the currents inferred, as 
above, from the calculated northward current den- 
sity. A more thorough evaluation of the field- 
aligned currents will be presented later. 

Acknowledgment . The Editor thanks J. R, 
Burrows for his assistance in evaluating this 
brief report. 
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Abstract. The ioii drift technique of measuring ionospheric electric fields is compared to two other 
simultaneous measurements. Rocket measurements in the evening auroral oval are checked against a 
dual probe, also in situ, and the ground based STARE auroral radar. The technique is explained thoroughly 
as well as tested for its dependence on mass. 

Two evening auroral oval conditions were observed from Andoya, Norway in January and February 
1977. The first flight, 18:1005, measured electric fields over a quiet pre-midnight discrete arc. An interesting 
plasma convection reversal wus observed poleward of the arc. The subsequent flight, 18:1004, samples 
a breal up phase aurora nearer local midnight. 


1. Introduction 

The intent of this paper is to compare the ionospheric electric field inferred from 
rocket in situ ion drift measurements with those obtained from dual probe obser- 
vations made aboard the same rocket and simultaneous ground-based auroral 
radar soundings during two auroral (lights. Such a comparison is necessary to give 
confidence to the technique of obtaining ion drift, hence electric fields, from measure- 
ments of the ion distribution function in the ionosphere. Hanson and Heelis [1] 
have reported difficulty in making ion velocity measurements with "various kinds 
of sensors" flown on sounding rockets. Moreover, the first reported use of directional 
intensity ion measurements from 0 to 5 eV by Whalen et al, [2] to obtain drift velo- 
cities aboard a sounding rocket flown into an aurora gave unprecedented large 
transverse electric fields (100-200 mV/m) and strong flow along the magnetic field 
(~2 km/s). 

Multiple electrostatic analyzers detecting ions from 0 to 5 eV/charge are used to 
measure the thermal ion distribution in energy and space. The ion flow velocity 
relative to the sensors is obtained by fitting the data to a shifted (due to vehicle 
potential)-drifting Maxwellian. This technique requires knowledge of the ion mass. 
One of the analyzers flown on each of the two auroral flights to be discussed in 
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this paper was capable of measuring the relative ion abundance of masses 16 ± 2 
AMU and 30 ± 3 AMU. The sensitivity of the technique to the ion mass is assessed 
using these measurements. 

A crude test of the ion drift technique has already been reported by Morgan and 
Arnoldy [3] where data from one detector were used to obtain a distribution function 
over the duration of the Helm III rocket flight. The inferred electric held compared 
favorably with simultaneous Chaianika radar observations. This (light was not 
during an auroral display, therefore, the long period (200 seconds) of data collection 
to obtain the ion distribution was acceptable. However, since the density dependence 
on altitude as measured by the ground radar was input into the fit routine, the in 
situ observations were not totally independent. 

Finally, a comparison of the rocket dual probe measurenien.s and the ground 
auroral radar for one of the flights to be discussed has already been made [4], Very 
good general agreement was obtained when the electric field was greater than 20 
mV/m, the apparent threshold of the radar technique. A description of the dual 
probe experiment aboard the rockets and the auroral radar technique has been 
discussed in the report by Cahill el al. [4], hence will not be repeated here. 


2. Ion Drift Instrumentation 

The detection of the drift of ions was accomplished by multiple electrostatic analy/crs. 
These analy/ers had 90 cylindrical electrodes (R = 2 cm, AW « 2 mm) and employed 
Galileo Model 4025 Chaanellron electron multipliers as sensors. Ion spectra from 
0 -•> 5 eV in 32 energy steps were measured every 25.6 ms with an energy resolution 
of 8 The acceptance cone of the analyzers was 2x7 full angle, larger than the 
~ 1.5 angular rotation of the spinning rocket per energy step. Two analyzers were 
!lown, one viewing up the Held at 30 with respect to the spin axis of the rocket and 
the other at 90 . The analyzer design was similar to those flown aboard the Echo 
111 night [3j. 

A third analyzer was flown, viewing also at 90 with respect to the rocket spin 
axis, to measure ion masses or 16 ± 2 AMU and 30 ± 3 AMU, This instrument 
used the above described analyzer for energy selection with a low-resolution magnetic 
mass spcetomeler added between the electrostatic analysis and the channeltron 
detector [ 5j. The mass spectrometer consisted of a permanent magnet and post- 
accelerator. The 16 and 30 AMU mass separation was achieved by post-accelerating 
the 0 5 eV ions selected by the electrostatic analyzer by the proper amount (153 
volts for 16 AMU and 84 volts for 30 AMU) in order to be detected after 90 deflection 
in the magnetic field by the channeltron sensor. The post-accelerator was set to pass 
0 + (16 AMU) 9Q" 0 of the lime and 30 AMU masses (predominantly MO* at rocket 
altitudes of 200 250 km) the remaining time. The mass resolution of the instrument 
was 20" o< sufficient to separate the two mass groups. 

The geometric factor for the three detectors was obtained by laboratory calibration 
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using a simulated omnidirectional ion beam [5,6]. The electrostatic analyzer had 
a geometric factor of 9.6 x 10 f ’cnr sr eV eV and the mass analyzer detector was 
somewhat smaller at 4.7 x 10 ''em J sr cV eV. 


w. 

3. Data Reduction 


a 


a 


u 


The primary objective of the data reduction was to determine the bulk llow velocity 
vector. Additional results include ion temperature, ion density and rocket potential. 
Each data point was assigned a vector direction (0, tf>) in geomagnetic coordinates 
(?. along B) and an energy (0 5 eV). Data points were accumulated for two seconds 
to give adequate representation of the ion distribution. An iterative least-squares 
lilting routine was applied to this collection, resulting in the energy distribution 
function parameters of the observed data. Before discussing the fit routine further 
we would like first to discuss two presentations of the raw data which have proven 
valuable to help visualize some of the plasma parameters obtained by the lit analysis. 

The measured energy spectra are a function of six independent variables: energy, 
particle pitch angle and azimuth, and the three components of th» rocket ram velocity. 
The parameters needed to determine the electric fields are the magnitude and direc- 
tion of the E x B drift vector. The magnitude of this velocity (translated to energy) 
is nearly impossible to observe in the raw data since the dominant contribution 
to the energy peak of the ion spectra is due to the drop through the — 1.5 V rocket 
potential. The energy contribution due to the E x B drift is only about a quarter 
of an eV. The energy contribution due to the rocket ram velocity is also approximately 
n.quarter of an eV. The direction of relative velocity between the rocket and plasma 
(/: x B drift plus ram) can, however, be readily viewed in the raw data using a 
three-dimensional mapping program plotting detector counts vs. ion pitch angle 



tJf Pig. I. Example of data display (Plight 18:1005, time 275 seconds) of thermal ion llow velocity relative 
m to rocket, Axes arc pitch angle and azimuth, vertical is analyzer response in counts integrated over energy. 
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N'S. ion azimuth (Figure 1). In this figure the dimension front to rear is linear in 
pitch angle, O ' 180 , Azimuth is plotted from left to right, 0 360 , measured 

from north through east. The vertical axis is average counts integrated over energy. 
Essentially what is displayed is the velocity directions observable by the 30 and 
90 electrostatic analyzers. Thus the How direction can be visually located by a 
peak and compared with the answers obtained by the multiparameter Maxwellian 
fit to the data. The data shown in Figure 1 were taken on the downleg of flight 1 8 : 1005 
at a flight time of 275 seconds. As will be shown later, the least-squares fit predicted 
this time to have a generally southward electric field. This field results in an eastward 
E x B drift. Since the rocket ram is included in Figure 1 it must be subtracted. 
The ram at 275 seconds was south as the rocket, was first nearly straight north 
The rc&olt is an eastward flow combined with a comparably equal south ram velocity 
or a southeast peak (~ 135 ' azimuth). This can be seen in Figure 1. These plots arc 
therefore most useful for locating the azimuth of the flow. The limited pitch-angle 
coverage for two detectors is evident in the figure. 

The same plotting package was used to plot detector counts vs, energy vs. pitch 
angle (Figure 2). In this figure are plotted a series of energy spectra ranging in pitch 
angle. Consider the energy dimension as a radius arm originating at the center 
back. From this point along the back to the left is O' pitch angle. The radius arm 
pivots about its origin in a right-hand sense with increasing pitch angle. A radius 
arm starting at center back pointing right would be antiparallel to B or 180' pitch 
angle. The vertical axis is again counts averaged over azimuth. This plot shows the 
spectral shift in energy due to the vehicle potential. It also shows that the flow is 
mostly perpendicular to B, in which case the peak flux would occur at 90', with 
ram velocity contributing a component greater than 90'' pitch angle for this flight 
time of 275 seconds. Figure 2 is especially useful for detecting noise in the energy 


Flight 16 1005 

Pilch Angie vs Energy vs counts 



Fig. 2, Similar to Figure 1, showing pilch angle and energy in polar coordinates. A radius arm length 
starting at the rear center indicates 0-4,5 eV ion energy and its angle from the left rear corner (labelled 
0 ) is the pitch angle. Vertical is analyzer response in counts. The base line semicircle is 4,5 eV, 
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sweep. Figures like 1 and 2 have proved invaluable for initial analysis and parameter 
confidence. 

Progressing now to the fit routine, the original thermal ion velocity distribution 
is assumed to be Maxwellian: 


m = n 


InkTJ 




( 1 ) 


where the prime indicates that the independent variable is outside the rocket Debye 
sheath. The differential directional intensity J (E') is related to the velocity distri- 
bution according to: 


2 E' 

J(E') = /(f') 

nr 


particles 
s cm 2 sr eV 


( 2 ) 


In passing through the rocket sheath, the unperturbed Maxwellian is shifted according 
to Liouville’s theorem: 


J{E)/E = J{E')/E'. 


( 3 ) 


Finally, including a relative velocity between the rocket frame and the plasma, 
the final distribution recorded by the detector becomes: 


J(E) 


2 n 
in 2 


m 

2r.kT 




~([2(E - U)/m ] 1/ 2 f>- v 0 
2kT/m 



counts 

G 0 AtE 


( 4 ) 


where C 0 = geometry factor = 9.6 x JO" 6 (cm 2 ster eV/eV) energy analyzer or 
4.7 x I0“ 6 (cm 2 ster eV/eV) energy mass analyzer, E= selection energy (eV), 
At = accumulation time = 8 x 10“ 4 (s), n = density (cm -3 ), in — mass (g), T = 
temperature (K), U — energy due to rocket potential (eV), f = particle direction, 
v r — rocket ram velocity (cm/s), v 0 = drift velocity (cm/s). 

The data are fit to the above function by iteratively adjusting the six parameters 
(3 components of v 0 , T, U, and n) and finding the minimum of the least-squares 
function [7], Two seconds of data were used for each computation representing 
blocks of data taken over 78 spectral sweeps or ~2500 data points accumulated 
over 10 spins of the rocket. A typical pitch-angle vs. azimuth coverage by the two 
electrostatic analyzers per spin can be seen in Figure 1. Pitch angles greater than 
130" were not sampled. A T coning angle of the rocket slightly extended the data 
cos 3rage over pitch angle during the two seconds of accumulation. 

The mass of the ions detected is a parameter in Equation (4), The third detector, 
capable of measuring ion spectra for masses 16 mid 30 obtained the ratio of these 
masses at various altitudes. Figure 3 gives the measured ratio of the number of 30 
AMU ions to the total number of ions as a function of flight time and altitude for 
Flight 18:1005. The rocket was over an 8 kR (3914 A) arc from 140 to 155 seconds. 
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Fig. 3. Measured ratio of "0 AMU ions to total as a function of time (altitude). 


Beyond 155 seconds there was no further auroral luminosity, although electrons 
less than 300 eV with directional intensities exceeding I0 H (cars sr kcV)“ 1 were 
measured until 280 seconds flight time. The factor of two decrease and subsequent 
increase appears to be an altitude dependence and is roughly consistent with the 
altitude dependence of 30 ± 3 AMU tons given by Bauer [8]. 

Two different fitting routines were used with the data from Flight 18:1005; one 
assumed that the mass in Equation (4) was () , (1 and the other fit the data to a 
Maxwellian for 16 AMU mass plus a Maxwellian for 30 AMU mass using the 
measured composition ratio. Table 1 is a compilation of fit parameters using the 


TABLE) 

Fit parameters for two fitting routines using data from Flight IS: 1005. 
Method AM U 30 = 0.35 total ALL O * 


•on 

(m s) north 

-256 

±60 

-289 

± 64 

f-’SX 

for 

(m s) east 

825 

±34 

919 

± 34 


*OH 

(m s) along B 

-257 

± 66 

— 361 

±- 86 


Temp 

IK] 

1520 

+ 66 

1670 

± 72 


V 

(volts) 

1.72 

+ 0.004 

1.72 

± 0.006 


n 

, ill - 1 ) 

3.20 x 10 J 

± 800 

3.07 x I0 J 

±200 

f 
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two methods for one interval of data. These data were acquired at about 280 seconds 
flight time at an altitude of 215 km 

Generally there is about a 10 15"„ discrepancy between the two methods, which 
is comparable to the error in the parameters. The conclusion is that mass analysis 
is not essential to the determination of electric fields by this technique in that the 
accuracy of the spectral measurements cannot resolve the effect of an admixture 
of mass 16 and 30 ions. 








hig 4 DMSP photograph, courtesy of t red rick R'ch. Rex is College Research Center and Air Force 
Geophysics t ab taken at about eight minutes after I light IX Iihi5 Discrete trajectory is indicated 
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4. Oulu Presentation 

The object of both (lights was to vise the high spatial and temporal resolution of 
rocket techniques to study the premidnight auroral oval environment. In order to 
make the observation more complete, the two launches wore at different magnetic 
local limes, 2118 and 2250. 

Fluiln IS: 100? 

This llight was launched over a elassie premidnight discrete are. The are was stable 
and of wide longitudinal extent shown by the DMSP photograph of Figured taken 
approximately S minutes after the llight. The 'inverted V' in the measured electron 
spectrogram of Figure 5 locates the are between 135 155 seconds flight lime. The 
3914 A light emission by this are as measured by onboaixl photometers and from 
the ground was about 8 kR. The are was in the region of an enhancement of the 
eastward eleetrojet current, as calculated using conductivity derived from electron 
flux measurements and the measured electric field [4,9]. 

Figure 6 gives the simultaneous measurement oft' perpendicular to the magnetic 
field by the ion drill experiment and the dual probe experiment of the University 
of Minnesota. Figure 7 is the comparison of the dual probe with the STARE ground- 
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based radar electric field reported earlier [4]. General agreement can be seen between 
n the three methods, particularly in direction. There is a systematic discrepancy in 

il electric field magnitude between the ion drift method and both the dual probe and 

STARE measurements. This may be remedied with belter pitch angle coverage or 
r - the distribution function. (Statistical error in the velocity vector by the ion drift 

|| method was about 28°o. This would account for the general disagreement in magni- 

a tude between the dual probe and ion detector measurements of Figure 6.) Beyond 

290 seconds, the count rate of the ion detectors becomes very low due to the change 
in ram direction during the downleg of the trajectory. Since this was the time when 
" the auroral particle precipitation decreased (Figure 5) the ion density was reduced. 
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Fig, 7. From Cahill cl al. [4], comparisons of electric field measurements (Flight 18:1005) by the dual 
probe and the ground based STARE radar. The STARE E field vectors have solid points while the rocket 
data point in each frame has an open point. Flight time indicated is the time at the beginning of the STARE 

integration period. 


The STARE radar signals dropped to a low signal-to-noise ratio at £-region altitudes 
conjugate to the rocket at this time as well. 

Geophysically, a significant event in these data was the (low reversal at about 
200 seconds. The electric field reversed from generally poleward across the arc and 
up to the boundary at 200 seconds to equatorward after 200 seconds. The poleward 
field was about 25 mV/m while north of the boundary a larger southward field was 
observed. Note that the STARE observations indicate the reversal to be a spatial 
structure. These observations are consistent with the transition from sunward to 
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antisumvard plasma convection measured by polar-orbiting satellites [10,11] at 
the edge of the polar cap, 

Flight IS: 1004 

This flight was launched closer to midnight and duting the breakup phase of aurora. 
Although there was strong westward electrojet activity during the flight, the auroral 
display had faded to very low levels as the rocket was near the top of its trajectory. 
As might be expected, the aurora became very brilliant after the rocket flight, reaching 
intensities as high as 50 kR in 5577 A. The onboard 5577 A photometer recorded 
patches of aurora not exceeding 2 kR during the flight. 



Fig. 8. Comparisons for Flight 18: 1004 of the electric field vector as measured by the ion drift detector 

and the dual probe instrument. 




394 


I . J. /AN! TTI IT At . 


The comparison of electric field inferred from the ion detectors with the dual- 
probe results is given in Figure 8. Again, the magnitude of the field from the ion-drift 
detector is larger than that obtained by the dual probes but the directions are in 
good agreement near the middle of the (light. There is a systematic discrepancy in 
direction between the ion drill and dual probe measurements at the end of the 
(light. The STAKE measurements during the flight alon» with the dual-probe data 
are given in Figure 9. At this time, the rocket trajectory and the STARE radar were 
not eo-loeated. The westward boundary of the radar field of view is indicated in 
Figure 9, The measurements to the east of the rocket trajectory support the in- 
flight electric field measurements, 


o 

(S> 


O) 

XD 

ZD 


Flight 18:1004 

Store E-Fields vs. Rocket E* Fields o-» 

50 millivolts /m * 

Integration Time (sec): 20 



Pig. 9. STARE i; Reid observations in comparison to dual probe measurements, Flight 18:1004. The 
westward boundary of the radar field of view is indicated 
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There is tin interesting dynamic pattern that seems to be drilling across this 
section of the oval, The sequence is a slight northward turning of the westward 
electric field followed by a transition southward. This disturbance appears to be 
drilling eastward as it is first observed by the rocket and later seen in the STARE 
field of view. Note in Figure 9 how the rocket measurements at 197 seconds flight 
time are seen by the STARE at 217 seconds. 217 seconds rocket observation viewed 
at 237 seconds by STARE, and so on, A rough calculation indicates the eastward 
flowing perturbation to have a drill of slightly greater than 1 km s. This is a commonly 
observed electron drift velocity in the westward clcetrojet region although the 
predominantly westward electric field observed at this lime indicates a southward 
drill. 

In conclusion, the ion-drift technique of inferring ionospheric electric fields from 
a rocket platform appears to give results consistent with the more accepted dual 
probe and auroral radar observations, The large error that we have assigned to the 
ion-drifl technique for the measurements presented in this paper is largely due to 
the relatively low count rates obtained. The instrument was designed to operate in 
ambient densities of 10* to 10 (> cm" 3 . The densities encountered in the two flights 
were about 10 4 cm 3 as obtained from the fit routine. These densities arc quite a 
bit lower than those normally observed, especially in the auroral regions. This 
problem is still being investigated. One should note, however, that the drift velocity 
vector is determined solely from position in 4rt stcradians for direction and the 
location in energy of the peak II ux for magnitude. Thus, the drill velocity is inde- 
pendent of density. Recent calibration tests have indicated dillerent efficiencies for 
low energy electrons or low energy ions impinging on Channeltron cone surfaces. 
This will effect G 0 in Equation (4) and, therefore, n, however nG 0 is unchanged. 

Finally, we would like to briefly describe the other parameters measured by the 
ion drill detectors, although there is no way to compare them with an independent 
measurement either from instrumentation aboard the rocket or on the ground. 
For Flight IS: 1005, the vehicle potential rose to -1,7 volts in the arc whereas 
outside the are it remained near — 1.5 volts. Inside the arc, the ion temperature was 
inferred to be ~ 3000 K. whereas outside, a value of 1500 K was obtained. For the 
fading aurora of Flight 18:1004, the temperature was about 2000 K. Near the field 
reversal seen in Flight 18:1005 at 200 seconds flight lime, the ion detectors sensed 
an enhanced flow of ions out of the ionosphere resulting in an upward field-aligned 
current of - 2 /(A m 2 . This upward current svas comparable to that carried by the 
auroral electrons producing the arc. This point will be further discussed in a sub- 
sequent paper when the arc of Flight 18:1005 will be modeled using all the field 
and particle measurements made aboard the rocket. 

Future flights involving the ion drift technique will have additional plasma sensors 
to measure electron density and temperature and vehicle potential so that a direct 
comparison can be made for these parameters as we have done here for the electric 
field. These flights will also use four ion drill detectors giving four simultaneous 
pitch angle measurements of the distribution function per azimuth increment. 
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The Hot Plasma Environment and Floating Potentials of an 
Electron-Beam-Emitting Rocket in the Ionosphere 
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Willi an extensive array of particle sensors the plasma environment surrounding the Echo III accelera- 
tor payload is studied. From measurements of the thermal ion spectrum, negative payload potentials ref- 
erenced to the unperturbed ionospheric plasma arc obtained. Multiple detectors determined the electron 
population from a fraction of an eV up to 40-keV energy. An intense electron population extending up to 
the energy of injected electrons by the accelerator is produced when the acceleration is turned on. The 
energetic tail of this population is returned to the payload principally from directions in which the beam 
was fired, with the most intense fluxes coming from the atmosphere for downward injections of the beam, 
The atmospheric-scattered beam and secondary electrons are called ‘quick echoes.’ Electrons of energy 
less than several keV down to the detector threshold (0,1 cV) are called the suprathermal component. 
These electrons are produced isotropically around the payload during gun firings and decay away in ap- 
proximately 32 ms. The largest directional intensities of this component arc observed at the higher alti- 
tudes. Quick echo electrons are also observed to produce suprathermal electrons when they encounter the 
payload. The mechanism by which the suprathermal electrons arc produced is discussed but remains un- 
known at the present time. The hot electrons surrounding the accelerator payload during gun injections 
bring sufficient charge to the payload to neutralize it provided the loss of charge by secondary production 
on the payload skis' is small, presumably owing to a positive payload floating potential during injection, 
Since the hot population exists for tens of milliseconds after the gun urn off, it results in driving the pay- 
load up to 4 volts negative during this time. Secondary production on the payload skin apparently pre- 
vents much larger negative potentials after gun turnofT when the payload is immersed in the hot popu- 
lation. Quick echo electrons creating suprathenrmi electrons around the payload also drive the pnyload 
to a few volts negative. This is an important consideration when discussing the origin of the suprathermal 
electrons, since electrical discharge and beam plasma discharge mechanisms do not apply to the quick 
echo beams. 


Introduction 

The principal purpose of this paper will be to analyze the 
particle distribution functions measured on board the Echo 
III electron beam experiment in order to deduce the floating 
potentials of the vehicle and to identify the source and nature 
of the hot plasma surrounding the rocket. The subject of pas- 
sive vehicle potentials in space has been much discussed [Al‘- 
peri et al„ 1965; Kasha, 1969], and it is well known that the ap- 
proximate floating potential of a sounding rocket in the 
ionosphere may be 1 V negative. But if the rocket is the source 
of an injected electron beam, such as the University of Minne- 
sota Echo series and numerous other artificial electron beam 
experiments, then the equilibrium becomes violently dis- 
turbed, and a net current must be drawn from the surrounding 
ionosphere equal to the electron beam experiments, then the 
equilibrium becomes violently disturbed, and a net current 
must be drawn from the surrounding ionosphere equal to the 
electron beam current carried out of the local region consid- 
ered. Generally, it is expected that during electron beam in- 
jection the vehicle potential will rise positively to maintain 
current balance, An important result of the beam injection is 
that the ambient plasma becomes heated and may be highly 
turbulent [Cartwright et al., 1978], and a new plasma environ- 
ment is in effect created by the beam itself. Even after the end 
of electron pu’vs the surrounding medium is found to remain 
appreciably heated for many tens of milliseconds. The float- 
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ing potential of the rocket may be considerably altered in this 
hot plasma medium and be much more negative than in the 
normal ionosphere. 

The term ‘vehicle floating potential’ will be defined opera- 
tionally here along the following lines for the Echo III experi- 
ment: The ambient ion distribution unperturbed by the pres- 
ence of the rocket payload is assumed to be Maxwellian. The 
vehicle potential with respect to this unperturbed plasma is 
determined for the Echo III flight by measuring the thermal 
ion spectra reaching the rocket sensors after penetrating the 
sheath across which the potential drop occurs surrounding the 
payload. As examples, theoretical ion spectra for floating po- 
tentials of +0.5, and —1.2, and —4 voh re given in Figure 1 
along with the unperturbed (0 V) Maxwellian. These spectra 
fora density 2 X 10 J cm 1 and a temperature of 1500°K, val- 
ues which correspond closely to the average ionospheric val- 
ues during the Echo III flight, arc expressed m terms of the 
differential directional intensity reaching the detector. The ef- 
fect of a relative velocity of ±10 5 cm/s between the plasma 
and the detector on the rocket, also close to actual values, is 
given for the -1.2- and — 4-V spectra. Negative floating poten- 
tials are readily identified by the shift in the entire Max- 
wellian ion spectrum toward higher energy with a sharp mini- 
mum energy cutoff resulting from acceleration of the ambient 
ions in passing through the sheath, Actual measured jon 
spectra, as we shall show, are consistent with those theoretical 
predictions for negative floating potentials, and thus there is a 
clear procedure "or evaluating such negative vehicle charging 
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Fig. 1. Theoretical ion spectra ( ,< 0 + ions) for an F region temperature of 1500° and a density of 2 X 10 5 cm -3 . Spectra 
are given for payload potentials of +0.5, 0, —1.2, and —4 V. The solid curves are for zero rocket ram velocity, and the 
dashed curves for ±1000 m/s ram velocities. 


by simply measuring the low-energy spectral cutoff. All neces- 
sary details of the ion measurements and analysis for the Echo 
III experiments have been published by Morgan and Arnoldy 
[1978]. This technique has been utilized with success by De- 
Forest [1972] and Whipple [1976] in revealing large negative 
excursions of up to 10 kV for synchronous orbit satellites 
when exposed to the natural energetic electron component 
during substorm events. Recently, a special synchronous orbit 
vehicle has been successfully launched to study such questions 
[Steve/if, 1979]. 

In contrast, without knowledge of the unperturbed ion 
plasma parameters and the nature of ion drift motions, posi- 
tive floating potentials cannot be uniquely determined from 
ion spectra measurements, It might be thought that electron 
observations should reveal positive potentials also by a mini- 
mom energy or cutoff point in analogy with the ion measure- 
ments. However, scattering and thermalizing, rocket shield- 
ing, secondary production by energetic beam electrons and 
other effects seriously obscure these measurements; con- 
sequently, it has been difficult to interpret the results of re- 
tarding potential analyzers when applied to electrons [ Winck- 
ler, 1976]. Recently, by a combination of on-board particle 
and remote electric field techniques using a mother-daughter 
payload, some large positive vehicle potentials have been 
identified with the Norwegian Polar V experiment [Jacobsen 
and Maynard, 1978]. 

Instrumentation 

The thermal ion and suprathermal electron distributions 
surrounding the Echo III rocket were measured by elec- 
trostatic analyzers aboard the payload. The ion electrostatic 
analyzers consisted of 90° cylindrical electrodes (mean radius 
of 2 cm) and a ‘channeltron’ multiplier biased to post- 
accelerate the ions by 1500 V as the sensor. Thirty-two-point 
spectra from 0 to 5 eV were measured every 25.6 ms by two 
detectors positioned parallel and perpendicular to the payload 
spin axis. A more complete description of these sensors has 
been given by Morgan and Arnoldy [1978], The electron detec- 


tors measured the directional intensity perpendicular to the 
rocket spin axis. Two cylindrical analyzers with different geo- 
metric factors and energy ranges were flown. Electrons from 
100 eV to 15 keV were measured in 128 windows every 200 ms 
by an analyzer with a geometric factor of 4.6 X 10~ 3 Sr cm 3 
keV (keV)" 1 . The so-called ‘return current analyzer’ measured 
the electron spectrum from 0.1 eV to 1.2 keV in 128 windows 
every 200 ms, Although three data points were taken below ? 
eV, these are to be regarded with caution, since they require 
very small potentials (<150 mV) across the analyzer. The en- 
ergy sweep of the detector was accomplished by a 200-V po- 
tential across the electrodes decaying to 0 V in 200 ms. Volt- 
ages of 200 mV and higher could be substantiated by 
measurements of the voltage decay curve during the flight; 
lower voltages remain uncertain. This detector had a very 
small geometric factor, 2.5 X 10' 6 sr cm 2 keV (keV) -1 , in order 
to sense the neutralizing current during gun firings. 

Finally, electrons from 10 to 40 keV were sensed by an ana- 
lyzer having eighth-spherical electrodes making a 256-point 
spectral measurement every 200 ms, A large geometric factor, 
2 X 10 -3 sr cm 2 keV (keV) -1 , was used to attempt a measure- 
ment of gun electrons echoing from the opposite hemisphere. 
The combined set of three electron detectors covered the en- 
ergy range 0.1 eV to 40 keV and made directional intensity 
measurements from 10“ to 10 13 electrons (cm 2 s sr keV)“', Fig- 
ure 2 is a schematic showing the acceptance directions of the 
ion and electron (auroral) detectors and the injection direction 
of the electron accelerator with respect to the local magnetic 
field. The Black Brant motor remained with the payload, and 
the entire vehicle had a spin period of 0,365 s and a coning pe- 
riod of 43 s. 

Measurements of Vehicle Floating Potential 

The electron gun flown aboard the Echo III flight was fired 
according to a preprogramed sequence of various electron 
energies (29-43 keV) and pulse duration (2, 8, and 32 ms), 
with periods of silence between gun sequences. The ion elec- 
trostatic analyzers were used during these undisturbed periods 
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Fig 2, Schematic showing the tiemr'taaee directions of the ton 
and electron (auroral) detectors and the ■ • etion direction of the elec- 
tron accelerator with respect to the local magnetic field, The coning 
half angle was 12°. 


to infer the F region electric field from the B x B convection 
drift imposed by the electric field upon the ambient ion popu- 
lation which appears ns a relative velocity between the detec- 
tors and the plasma in addition to the rocket ram effect. The 
ion detector measurements during and just following gun pul- 
ses, injecting about 70 mn of 40*keV beam electrons into the 
ionosphere, will be discussed in this paper. 

Figure ,1 gives three successive ion spectra measured by the 
detector looking parallel to the rocket spin axis. Energy in- 
creased with time during the spectral sweeps such that 1 eV *« 
5 ms. The first lin’d third spectra (sweep N nnd N 2) were 
measured betsveen gun firings nnd represent ambient con- 
ditions. An 8-ms duration electron injection (heavy hori- 
zontal line in Figure 3) occurred during the second spectrum 
(TV + 1 sweep) just at the energy (time) where the ambient 
population would normally have been detected. Some tons 
were detected during the gun pulse showing that the rocket 
potential even during beam injection occasionally was close to 
0 V. One should note that the voltage applied to the gun was a 
l-kllz. square wave resulting in a partial or complete beam 
dropout each millisecond. The ‘spikey’ appearance of ions 
during, gun-on periods ns seen in Figure 3, center panel, may 
reflect this fact. Several milliseconds after the gun turned oft', 
the mam ion population was found between 3 and 5 eV, Note 
that for the next spectral sweep, (.N + 2) after 25.6 m, the ion 
population had returned to its usual position. The nature of 
the observed spectral perturbation resulting from electron in- 
jection depended upon when the gun was fired in the spectral 
sweep. Figure 4 is an example of a gun firing just after the en- 
ergy analyzer sweep pasrztd through the region of the ambient 
population during sweep N. The ion population perturbed by 
the gun firing is now located between 4 and 5 eV nnd was 


measured in the same sweep after the gun turned oil’. In the 
next sweep, (A 4- 1), the ion spectrum had not yet quite re- 
turned to equilibrium conditions. 

If we interpret the leading edge and low-energy cutolf of 
the ion spectra as the vehicle floating potential, we see in the 
examples presented above that immediately following gun 
pulses this potential is observed to be several volts more nega- 
tive than its nominal -1,2-V value. Since the ion electrostatic 
analyzer sweep was randomly phased with respect to gun 
pulsing, it has been possible to construct a time profile of ion 
spectra cutolf nnd therefore vehicle potential before, during, 
and after gun pulsing using about thirty 2-ms gun pulses near 
rocket apogee. The data are shown in Figure 5 and have been 
separated according to whether the gun electrons were fired 
down info the atmosphere (« ** 70°) or away from the atmo- 
sphere (n « 110®). The payload tended to reach its largest 
negative potential ~ 5 ms after the end of the gun pulse and 
then returned to the equilibrium potential of —1.2 V by 10 ms 
for injection upward, For downward injection the negative 
payload potential did not decay back to the equilibrium value 
until about 20»ms after the gun firing. This long-enduring po- 
tential for downward injections is apparently due to the im- 
mersion of the payload in a population of backscattered beam 
and secondary electrons from the atmosphere reaching the 
rocket approximately 12 ms after gun firing. In a detailed 
study of particle dynamics in the vicinity of a beam-emitting 
rocket, Winckkr ft til. [1475] have referred to such particles as 
quick echos. Figure 6 gives an example of such a quick echo 
(heavy line) registered by one of the Echo 111 electron detec- 
tors following an 8-ms injection of 24 keV electrons down- 



ing. 3. Measured ion spectra for three consecutive energy Sweeps. 
The dark horizontal bar in sweep N + I corresponds to the time of a 
gun firing. The spectral sweep time was 25,6 ms, sampling from 0- to 
5-eV energy. The gun injected electrons at 47° pitch angle, and the 
ion sensors detected ions having a pitch angle of 53°, 
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Fig, 4. Ion spectra for two energy sweeps, Cun-on time is ngiiiii in- 
dicated by the dark horizontal bar. 


ward into the atmosphere, Tlte detector wits in the process of 
scanning in energy from 19 keV to 14 keV during the echo. 
The maximum response to the returning electron beam is de- 
layed 12 ms from the center of the gun pulse, which is approx- 
imately the time required for a 30-keV electron having a 70“ 
pitch angle at 270 km (which lies in the loss cone) to travel to 
the 100-ktn altitude level, undergo atmospheric scattering and 
magnetic reflection, and return to the rocket. The observed 
energy is approximately that expected for backscattered 29- 


keV beam electrons [ Winckler et al„ 1975). The light line in 
Figure 6 shows suprathermal secondaries of 30=150 eV near 
the rocket generated both by the injected beam and by the 
quick echo electrons. 

One can interpret the increase in negative payload potential 
following gun pulses shown in Figure 5 as due to a rise in the 
temperature of the r t Minding electron population which we 
shall now proceed * tbs \ss. Langmuir probe measurements 
made aboard the Echo t. night {Cartwright ft «/., 1978) have 
shown that the electron accelerator heated the electrons of the 
background plasma to IQ,000°K, Cartwright ft al. (1978) sug- 
gested that the high-temperature component of the back- 
ground plasma represented secondaries produced by the ac- 
celerator beam. 

Tm- Hot Electron Population 

The Echo 111 data are unique among such experiments in 
that the electron spectrum reaching the payload was measured 
over a very large energy range, between 0.1 eV and 40 keV. 
The hot electron population surrounding the payload, neces- 
sary information for payload electrodynamic studies during 
electron injection, is therefore known. Spectra measured at 
various times with respect to 8-ms duration gun firings and a 
background spectrum due to natural electrons are given in 
Figure 7. This figure gives the dill’erential directional intensity 
of upward moving electrons generated by downward electron 
beam injections of 8-ms duration, at nominally 40-keV energy 
and 70-ma beam current. Since the detectors required 200-ms 
to make a complete spectral measurement, data were com- 
piled over a (light time of ~ 150 s (fifteen 8-ms gun pulses) in 
order to obtain the various spectra of Figure 7 measured at 
the indicated times with respect to gun firings, Limited energy 
coverage at the desired time with respect to gun pulsing ab- 
breviated some of the spectra at low and high energies. Ex- 
amples of actual data points are given for the upper and lower 
curves but not for the remaining spectra to avoid confusion in 
the presentation. 



Fig, 5. Floating potentials of the Echo 1U rocket in the quiescent ionosphere at 250 km before beam injection (f < 0), 
during 2-nis beam pulses (quite likely fluctuating and generally unknown), and after injection, deduced from the low-en- 
ergy cutoff of the spectra of thermal ions, A*. 5 ms the rocket reaches an extreme negative potential in the heated plasma 
created by the beam. At IS ms another negative excursion is associated with quick echoes returned from below the rocket 
(only seen for downward injections). 
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Pig, 6, Directional electron intensities measured by two detectors 
during and following an injection of 29-keV electrons into the atmo- 
sphere. The energies sampled by the two detectors ns a function of 
time are given at the bottom of the figure. The heavy curve is a quick 
echo. The gun injection and the quick echo electrons both produce su- 
prnthermal electrons (light curvet. 

The low-energy portion of the electron spectrum maximizes 
during the gun pulse and then decays to background after in- 
jection. In contrast to this a significant increase in the flux of 
high-energy electrons occurs after the gun turns olf, reaching 
a maximum intensity 8 ms after the pulse has terminated. This 
flattening of the electron spectrum is due to the quick echo 
electrons discussed above filling in the region above 10 keV, 
Purely magnetically mirroring gun electrons are not a part of 
the quick echo beam because they cannot reach the particle 
detectors as the payload moves off the field line of injection 
[ Winckler et al„ 1975], Only those gun electrons which have 
multipk-Coulomb scattered in the atmosphere below the 
rocket to field lines and directions bringing them back to the 
payload will be detected in the quick echo; hence one should 
not be surprised at the ubsencc of a peak in the backscattercd. 
spectrum at the gun energy. The hypothesis that the spectral 
flattening is due to atmospheric bacicscattcred electrons is sup- 
ported by measurements made when the gun fired away from 
the atmosphere, where the spectral flattening docs not occur. 

The response of the detector which sampled energies be- 
tween 150 eV and 10 eV plotted in Figure 6 is a specific ex- 
ample of the time profile of the supralhermal population fol- 
lowing a gun injection. The data of Figure 6 are for an 8-ms, 
29-kV downward electron injection. For injections upward 
where no quick echos would be produced, the hot population 
in this energy range would decay much more rapidly than the 
~20 ms given in Figure 6. The secondary peak in the hot pop- 
ulation at 290.027 s in Figure 6 is apparently produced by the 


quick echo electrons. The supralhermal electrons in the sec- 
ondary peak in turn generate the negative payload floating 
potential maximizing at 15 ms after gun turnofl' ns given in 
Figure 5, while the more intense fluxes measured during and 
immediately following gun injection are responsible for the 
prompt (~5 ms) negative payload potential, 

Production of Negative floatino potential 

At equilibrium, any object in a plasma will adjust its float- 
ing potential as defined earlier to keep the total current flow- 
ing to the object equal to zero. The thermal ion measurements 
made between gun pulses aboard Echo 111 show that to ac- 
complish this, the payload floated at -1,2 V, accelerating ions 
to the payload and retarding the electrons in compensation 
for the diflercncc in the thermal ion and electron speeds. 
However, when the payload was embedded in the very hot 
electron population Following gun injections, we saw above 
how the floating potential exceeded 4 V negative. Since rather 
complete thermal and suprnthermal electron observations were 
made above Echo III, we can quantitatively discuss the am- 
bient — 1.2-V and the -4-V floating potentials during injection 
in terms of this supralhermal population. It should be noted 
here that the Echo 111 flight was made at night, so photoelec- 
tron effects can be ignored. 

Ambient Conditions, —1,2- V Floating Potential 

A multiparameter fit to the Echo 111 ion data [Morgan and 
Arnoldy, 1978] inferred an in situ ion temperature of 1500°K, 



Fig. 7. Measured electron spectra at vurious times with respect to 
the 8-ms gun injections. The spectra are compiled over a (tight time of 
150 s, or fifteen 8-ms gun injections. The spectra are for gun injection 
downward into the atmosphere and detection of upward moving elec- 
trons. Data points arc given for only two spectra to eliminate con- 
fusion in the presentation. 
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A density of 2 X 10 5 cm ' 3 at rocket apogee was measured dur- 
ing the flight by the Chatanika radar. Figure 1 gives the theo- 
retical, zero plasma potential Maxwellian for these values 
along with the shifted Maxwellian distribution as it would be 
observed on a vehicle surface whose floating potential was 
-1.2 V, which corresponds to the measured value. The shifted 
Maxwellian is obtained from the Q-V distribution by con- 
serving JjE according to Liouville’s theorem, where J is the 
differential directional intensity, having units of particles (cm 2 
s sr keV) Integrating the N sweep of Figure 3 over energy 
gives 1.5 X 10 10 ions (cm 2 s sr) \ which compares favorably 
with the corresponding value of 2 X 10'° ions (cm 2 s sr)"’ ob- 
tained by integrating the theoretical curve of Figure I. 

The thermal ion current density incident on the rocket skin 
can be estimated from the integrated directional intensity by 
multiplying it by the electronic charge and effective solid 
angle J2 over which the ions reach the payload. We assume 
that the rocket is surrounded by a sheath of 1.2-V potential 
difference. A thermal ,6 0'' ion moving approximately parallel 
to this sheath at an energy ofkT** 0.13 cV will be accelerated 
as it enters the sheath toward the rocket surface, gaining 1,2 
V. The maximum angle measured with respect to the detector 
view direction (normal to the sheath) that these ions will have 
as they reach the detector will be 18°. The effective solid angle 
52 is therefore 0.3 sr, and the ion current density becomes 
« 10 /ia m 2 . 

A rocket payload moving relative to the thermal ion plasmn 
will record greatly distorted ion energy distributions, as seen 
in the example in Figure I, for overtaking and receding veloc- 
ities of 10' cm/s. From an analysis of the distorted equilib- 
rium ion spectra the E x B drift has been inferred for the Echo 
111 data [Morgan and Arnoldy, 1978). However, the ion mea- 
surements presented here made by the detector viewing along 
the rocket spin axis experienced the least relative velocity near 
apogee, hence the approximate the zero-velocity case. The 
current given above is therefore the best estimate of the aver- 
age ion current density that can be made. A similar ion cur- 
rent is obtained by computing the charge swept up per second 
by the transverse rocket motion through the plasma, since the 
transverse rocket ram speed is comparable to the ion thermal 
speed. 

At equilibrium the ion current density should be bulanced 
by the thermal electron cirrcnt to the -1.2-V payload: 


nq 


X 10 4 



~WkT. 


0. 1 pa m 


-2 


for a density of 2 X 10* cm 1 and 1500°K temperature, This 
current is extremely sensitive to the temperature. The temper- 
ature need only be increased to 2500°K to produce a current 
equal to the ion current. In fact, a 2500°K electron temper- 
ature is reasonable since there was a diffuse aurozal region 
covering the eastward flight trajectory. Furthermore, the elec- 
tron background measurement given in Figure 7 shows a sub- 
stantial flux ofhot electrons above 100 eV, and thus the condi- 
tion = /electron seems attainable. 


Postbeam Injection Conditions, —4-V Potential 

The measured ion spectrum reaching the payload at -4-V 
potential (Figure 3, A r + 1 sweep) has quite a different shape 
than the theoretical shifted Maxwellian of Figure 1. However, 
the directional intensities integrated over energy in both cases 
agree to within a factor of 2. The measured value is 3 X 10 10 


ions (cm 2 s sr)" 1 compared with 6 X 10 ,n ions (cm 2 s sr)” 1 from 
Figure 1, The measured spectrum of Figure 3 is perhaps dis- 
torted because the detector energy sweep was not fast in com- 
parison to the time constant of the floating potential; hence 
the payload potential could have changed significantly be- 
tween energy steps, 

As discussed above, the negative potential will collimate the 
ions reaching the payload, and for -4-V the effective soiid 
angle is 0.10 sr. The ion current density at the payload be- 
comes 10 pa m” 2 for the theoretical directional intensity. Thus 
despite a higher directional intensity than in the -1.2-V case 
the collimation reduces the effective solid angle to give equal 
current densities. In fact, the ~IQ pa in 2 current density for 
both negative potentials is equal, as it should be, to the total 
thermal ion current density incident on the sheath around the 
vehicle given by 


nq X 10 4 


8 kT, 


mn i 


1/2 


** 1 1 pa m ’ 


for a density of 2 X I0 5 cm 3 and 1500°K temperature. For a 
thermal sheath interface area not too different from the rocket 
area (the Debye length at 260 km for ambient conditions is ~1 
cm) the above equality simply means that all the ions entering 
the sheath are attracted to the payload and collected by it, 
Integrating the hoi electron spectrum measured 8-ms after 
the gun pulse in Figure 7 gives the directional intensity of the 
electrons that reach the payload for the -4-V floating poten- 
tial, This integration gives 3 X 10" electrons (cm 2 s sr) ', re- 
sulting in an electron current density reaching the payload of 
1,5 rna m'' 2 , over 2 orders of magnitude greater than the ion 
current. The electron current can be made equal to the ion 
current if secondary emission produced by the hot electron 
bombardment of the payload results in only 1% net electron 
charge collection. We have further evidence for this from the 
Echo I rocket flight, During this flight a large aluminized my- 
lar disc was deployed to aid neutralization of the payload. 
During quick echos this collector was driven positive appar- 
ently as a result of secondary emission when bombarded by 
the quick echo beam and the hot population produced by it 
[ Hendrickson et ah, 1975). From the Echo III data presented 
in Figure 5 we know that quick echoes produced a hot popu- 
lation similar to that resulting from injections of the primary 
beam. One can therefore deduce that in the absence of sec- 
ondary electron emission from the puyload the hot electron 
distribution of Figure 7 would necessarily have driven the 
payload much more negative than the —4 V observed in order 
to repel the low-energy electrons and maintain current bal- 
ance at equilibrium. 


Payload Potential and Neutralization 
During Gun injections 

Integration over energy of the ‘middle of 8-ms gun pulse’ 
curve of Figure 7 gives 2 X 10 12 electrons (cm 2 s sr)" 1 or a cur- 
rent density of 10 ma m” 2 arriving at the payload during in- 
jection. If this current density was uniform over the entire 
payload area of 10 m 2 , a total negative current of 100 ma was 
collected, which agrees well with the 70-rna beam of electrons 
leaving the payload. 

For an ionospheric density of 2 X 10 5 cm' 3 and a temper- 
ature of 1500°K the thermal electron current density would be 
2 ina/m 2 , and the thermal current to the payload becomes 20 
ma. To neutralize the payload during a 70-ma electron injee- 
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Fig, 8. Segments of measured electron spectra during 32r tns gun 
tilings, Continuous lines are data from the same spectral sweep. At 
low energies the curves separate according to whether the gun is in- 
jecting electrons into the atmosphere or away from it. 


tion, a current amplification of 3,5 was therefore needed, and 
in fact one may consider the suprathermal spectra of Figure 7 
to be showing this amplification effect, Secondary emission 
would not be a factor in this case, since presumably the pay- 
load is driven positive during gun firings, thereby returning 
any secondaries produced back to the payload. It is important 
to bear in mind that the electron detectors measured particles 
having pitch angles between 40° and 140°. This combined 
with the fact that a large percentage of the collecting area 
could be reached only by (low transverse to B means that if 
the measured electrons neutralized the Echo III rocket, as sug- 
gested above, it was accomplished by cross B particle motion. 
No measurements were made to estimate the neutralization 
current parallel to B. Measurements made aboard the Echo 
IV flight by Israelson and Winckler [1979] using metallic col- 
lectors, however, showed that there was no appreciable differ- 
ence in return current collected along the field from that de- 
tected perpendicular to the field. 

The thermal ion observations presented above clearly show 
that a few milliseconds after the gun was turned oil', the pay- 
load potential recovered and in fact swung negative. There re- 
mains the question, however, as to whether or how much the 
Echo III payload charged positively during gun injections to 
attract the neutralization current. As mentioned before, the 
thermal ion data during gun pulses cannot measure positive 
potentials but do not rule out the possibility, The small geo- 
metric factor return current electron electrostatic analyzer was 
capable of measuring the energy of electrons returning to the 
rocket during gun firings. Portions of a number of electron 
spectra taken during several 32-ms gun pulses are plotted in 
Figure 8. Each line segment represents data taken during the 
same spectral sweep and for one gun pulse. The lack of any 
spectral cutoff in the curves of Figure 8 essentially leaves the 
payload potential question unanswered in that it might be in- 


terpreted as evidence for (1) a positive potential in which the 
electrons have undergone complete thermnlizntion before 
being detected, (2) a positive potential in which there was a 
strong electron source in the sheath (e.g., beam ionization of 
neutral gas surrounding the rocket, secondary emission from 
payload skin, or a plasma discharge around the payload), (3) 
no potential drop at all, or (4) a rapidly fluctuating potential 
not capable of being resolved by the measurements made 
every 3.2 ms. The tact that there is a discernible difference on 
whether the gun fires up or down indicates that without a 
thorough understanding of the source of the hot electrons sur- 
rounding an accelerator payload they are not a reliable mea- 
sure of vehicle potential. Ambiguity between two different 
techniques of determining vehicle potential aboard the Araks 
rocket utilizing electron measurements [Grinyauz and Schulte, 
1976] further substantiates this conclusion. Measurements of 
electric fields and electrons on a payload separated from the 
beam-emitting rocket have produced evidence, however, for 
large positive potentials in a region of space including the 
rocket during electron injections {Jacobsen and Maynard, 
1978]. 

Production of Hot Electrons 

The hot electron spectrum reaching the payload during and 
following electron injections undoubtedly has several com- 
ponents. One comnoncnt might be electrons accelerated to the 
payload during injections which may be thermalized by local 
collisions and the generation of secondaries. Another part 
could be composed of electrons scattered from the beam and 
the atmosphere (quick echos) and all the secondaries pro- 
duced by beam and scattered beam electrons, The Norwegian 
Polar 5 experiment clearly showed the production of a halo of 
secondary electrons around the beam by detectors on the 
mother payload physically removed from the daughter puy- 
lond carrying the electron accelerator [Maehhun et al., 1978]. 
To explain the increase in secondary electron fluxes with alti- 
tude as observed in the Polar 5 experiment, Maehlum et al, 
[1978] have suggested wave-particle interactions us a source of 
hot electrons at high altitudes. 

The spectra of Figure 7 are only for electron injection down 
into the atmosphere and electrons detected moving up away 
from the atmosphere. For other injection and view directions 
they are considerably different. In an attempt to unravel the 
dependence of detector response on energy, gun pitch angle, 
and detector pitch angle we have used the phase map presen- 
tation developed in the analysis of the Echo 1! data [ Winckler 
et al., 1975]. In this presentation the measured directional in- 
tensity is plotted on a pitch angle of gun injection versus pitch 
angle of detection grid. The magnitude of the directional in- 
tensity is four-level-coded by the size of the data point. The 
phase map of Figure 7 is made for the average directional in- 
tensity of electrons having any energy between 5 and 25 keV 
sampled during four time intervals: during any gun pulse, 0-8 
ms after, 8-16 ms after, and 15-24 ms after the pulse, Two of 
the three electron detectors flown contributed to this map. 
Since the detector energy sweeps were not synchronized to the 
gun firings, the energies sampled for nny gun pulse were ran- 
dom, Although electrons from 5 to 25 keV were accepted for 
the plot, for any given time segment the energy changed by 
less than 1 keV over the averaging time. Finally, injection at 
all altitudes were used in constructing the phase map of Fig- 
ure 9, 

The Figure 9 phase map substantiates the conclusion 
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! f ig. 9. Phase map (explained in text) for electrons between 5 and 
25 kcV, The size of the dot increases with the directional intensity 
measured. 


reached in the earlier analysis of the Echo II data | Winckkr et 
al.. 1975] that energetic electrons whose trajectories clear the 
payload return from directions in which the electron gun has 
just injected beam electrons. Quick echoes predominate and 
lie in the “inject down -detect moving up* quadrnnt of the 
phase map. A significant but still unexplained return of ener- 
getic electrons occurs when the gun injects up and the detec- 
tors view electrons moving down the magnetic field. 

The phase maps of Figure 10 are for the 0.1- to 1.2-kcV su- 
pralhcrmal component. The data for these maps are further 
selected according to time segments measured with respect to 
beam injection and the altitude at the time of injection. The 
two maps on the left are for altitudes below 200 km, and those 
on the right between 200 and 260 km. The top panels of each 
set include the time segment during the injection and the first 
8-ms segment after gun turnoff, The bottom panels are for the 
8- to 16- and 16- to 24-ms segments after gun turnoff. Average 
directional intensities during each time segments are plotted as 
in the phase map of Figure 9, 

The range of average directional intensities included on a 
given phase map is given in the Sower left corner. The maxi- 
mum intensities measured have an altitude dependence, in- 
creasing with altitude. The phase maps also show that as a 
function of time after gun injection the suprathermal electrons 
decay in intensity as the spectra of Figure 7 show. During and 
immediately following injections the suprathermal electrons 
show little dependence on injection and detection pitch angles 
as given by the top two phase maps, Apparently, the supra- 
thcrmal electrons form an isotropic cloud around the payload 
during this time. As the intensify decreases with time, how- 
ever, there is an indication in the high-altitude data of a 
longer time constant for the decay in the quick echo quadrant 
apparently due to their delayed production by the energetic 
quick echo electrons as seen in Figure 6. 

The altitude dependence of the suprathermal electrons first 
seen in the Polar V data [Mae Id am et al., 1978] has led these 


investigators to rule out secondary production by beam ioni- 
zation as a principal source of the suprathermal electrons at 
high altitudes. A beam-plasma interaction mechanism was 
postulated, since it would become more efficient at high alti- 
tudes and explain the observed altitude dependence. One such 
mechanism, beam-plasma discharge originally discussed by 
Getty and Smullin f 1963), has recently been observed in 
plasma tank experiments by Bernstein et al. [1979], The results 
of the plasma tank experiments, however, suggest that the 
Echo III electron gun was not capable of igniting the dis- 
charge [ Winckkr, 1979]; hence the interaction was probably 
not the source of the hot population, Moreover, the fact that 
quick echo beams (orders of magnitude less intense than the 
injected beam) also produced suprathermal electrons clearly 
requires alternative sources. 

Can one rule out the possibility of secondary production by 
beam ionization of the neutral gas around the rocket? The 
study of auroral light produced by the Echo HI beam as a 
function of altitude by Israelson and Winckkr [1975] shows 
that the neutral density surrounding the payload is not signifi- 
cantly enhanced above ambient values as a result of payload 
outgassing; hence the ionization must occur only in the am- 
bient neutral atmosphere. Since the neutral atmospheric den- 
sity decreases by 2 orders of magnitude over the altitude range 
130 -260 km, the ionization rate near the rocket will do like- 
wise. If equilibrium conditions arc assumed, then the decrease 
in the recombination coefficient as a function of altitude docs 
compensate for the drop in production rate to give nearly con- 
stant secondary densities ns a function of altitude. It is not nt 
all evident, however, that the equilibrium conditions should 
prevail during and immediately following gun injections. 

The observations of a visible light glow around the Zarnitza 
payload [ Lyachov and Managadze, 1977] suggest that an elec- 
trical discharged occurred owing to the high positive payload 
potential during injection. Such a discharge might produce 
the suprathermal secondaries observed by the Echo III elec- 
tron detectors during and immediately following injection. 
Photometers fiown aboard the Echo III flight, however, give 
no evidence for any light produced by such a discharge [Aj- 
raelson and Winckkr, 1975]. Moreover, this mechanism could 
not produce the suprathermal electrons associated with quick 
echos, since the payload potential was measured to be only a 
few volts negative at these times. Finally, the isotropic distri- 
bution of suprathermal electrons during injections requires a 
rapid thermalization of the hot electrons. Since strong wave 
emissions have been measured during gun injections [ Cart - 
wright and Kellogg, 1974], it is quite likely that wave-particle 
interactions must enter the picture to accomplish this and 
could also be instrumental in producing the suprathermal 
component. 


Summary 

1. From measurements of the thermal ion spectrum, nega- 
tive excursions of the payload potential referenced to the un- 
perturbed plasma up to -4 V were observed several millise- 
conds -after gun firing, 

2. A hot electron population existed around the Echo III 
payload during beam injections, decaying to background in- 
tensities in a few tens of milliseconds. The hot population ex- 
tended in energy up to the injected electron energy. The ener- 
getic electrons of this population returned to the payload 
principally from directions in which the beam was fired. From 
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Pig. 10. Phase map of suprathcrinal (0.1— 5.2 keV) electrons separated according to altitude (left maps are Tor altitudes 
below 200 km, and right maps arc for altitudes between 200 and 260 km) and time after gun injections. Top two maps arc 
for the time segment during injection and the first 8-ms segment after injection. The bottom maps are for the 8- lo 16-ms 
and 16- to 24-fttS segments after gun turnolf. The range of directional intensities plotted in each map is given in the lower 
left comer beneath each map. 


the atmosphere the energetic electrons were beam and second- 
ary electrons Coulomb-scattered onto the rocket field lines 
(quick echos) and were delayed in arriving at the payload af- 
ter beam injection by their transit time to go down to the at- 
mosphere and come back to the payload. The downward re- 
turn of energetic electrons for upward injections resembles the 
Echo 11 result and is not as yet understood. 

3. The suprathermai component (less than a few keV) of 
the hot population appears to have been produced isotropic- 
ally about the payload during gun injection. The intensity of 
this component increased with increasing altitude. 

4. The suprathermai component of the hot population was 
responsible for driving the payload negative several millise- 
conds after gun turnoff. 

5. Energetic quick echo electrons also produced a supra- 
thermai component when they intercepted the payload, result- 
ing in a second negative excursion in the paylond potential for 
downward injections. 

6. Setting the electron current collected by the payload 
equal to the measured ion current several milliseconds after 
the end of gun injections requires that the measured supra- 
thermai electrons were only 1% efficient in bringing negative 


charge to the payload. One might interpret this in terms of a 
rather large secondary emission from the payload when bom- 
barded by the suprathermai electrons, 

7, The Echo III particle measurements give no informa- 
tion about positive payload potentials during injections. If the 
payload were sufficiently positive to prevent the loss of charge 
by secondary emission, the measured hot electron intensities 
surrounding the payload were more than adequate to provide 
the beam neutralization current during injections, 

8. The positive dependence of the suprathermai com- 
ponent on altitude makes it difficult to interpret its origin in 
terms of beam ionization of the neutral atmosphere. Wave- 
particle interactions, beam-plasma discharge, and electrical 
discharge of a highly charged payload should be studied fur- 
ther. However, the suprathermai electrons produced by quick 
echo electrons could not have been generated by the beam- 
plasma discharge of an electrical discharge around the pay- 
load. Upon completion of this work we received a manuscript 
accepted for publication authored by Johnstone and Sojka 
[1980] dealing with the heating of ambient thermal electrons 
to temperatures as high as 10 5 “K. by the passage of a field- 
aligned beam of suprathermai auroral electrons. Since the 
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measurements were made at altitudes above 600 km, where 
collisional effects arc negligible, Johnstone and Sojka suggest 
that they represent evidence for beam-plasma interaction. 
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